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CHAPTER  I 


INTRODUCTION 

An  adaptive  array  is  an  antenna  that  can  change  its  pattern  in 
response  to  the  signals  it  receives  [1,2].  By  forming  pattern  nulls  in 
the  directions  of  undesired  signals,  an  adaptive  array  can  protect  a 
communi cation  system  from  interference.  The  Applebaum  adaptive  array 
[1]  uses  a  steering  vector  to  control  the  direction  of  its  main  beam. 
This  array  nulls  signals  that  arrive  from  angles  other  than  the 
direction  specified  by  this  vector.  The  LMS  array  [2]  uses  a  reference 
signal  to  control  the  main  beam  direction.  Both  Applebaum  and  LMS 
arrays  can  be  used  in  communication  systems.  When  the  desired  signal 
arrival  angle  is  known,  the  Applebaum  array  is  the  best  choice.  The  LMS 
array  is  appropriate  if  we  wish  to  track  a  desired  signal  whose  arrival 
angle  is  unknown. 

The  LMS  array  requires  a  reference  signal  that  is  correlated  with 
the  desired  signal  and  uncorrelated  with  any  interfering  signals  [31. 
Methods  for  deriving  suitable  reference  signals  have  been  developed  for 
several  modulation  types  including  binary  phase-shift-keying  (BPSK) 

[31,  quadrature  phase-shift-keying  (QPSK)  [4],  conventional  amplitude 
modulation  (AM)  [5],  and  frequency-shi ft -keying  (FSK)  [6,7], 

Experimental  systems  have  been  built  that  verify  the  feasibility  of 
each  of  these  systems  [3,4, 5, 7]. 
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The  output  of  the  LMS  array  is  a  weighted  sum  of  the  signals  at  its 
inputs.  The  LMS  weights  maximize  the  signal  to  interference-plus-noise 
ratio  (SINR )  at  the  array  output  [8],  Therefore,  in  much  of  the  early 
work,  the  output  SINR  was  used  as  a  figure  of  merit  for  adaptive  array 
performance  [8-13].  However,  a  few  recent  studies  have  used  the 
bit-error  probability  at  the  communication  system  output  as  a  measure  of 
performance.  Al-Ruwais,  Acar,  and  Compton  [12,13,14]  estimated  P(e)  at 
the  output  of  differential  phase-shi ft -keyed  (DPSK)  communication 
systems  that  included  LMS  arrays.  Their  estimates  were  based  on  the 
assumption  that  the  noise  and  interference  at  the  array  output 
contributed  equally  to  the  degradation  in  detector  performance. 
Furthermore  they  assumed  that  the  detector  performance  was  not  dependent 
on  the  spectral  distribution  of  the  noise  or  interference. 

In  this  report  we  examine  the  performance  of  digital  communication 
systems  that  use  LMS  adaptive  arrays.  We  consider  the  performance  with 
BPSK,  QPSK,  binary  DPSK,  and  binary  FSK  modulation  of  the  desired 
signal.  We  examine  the  effects  of  single  tone  (i.e.,  CW)  interference 
and  bandlimited  Gaussian  noise  interference  on  these  systems.  We 
identify  key  operating  parameters  and  calculate  their  effects  on  system 
performance.  These  parameters  include  the  received  signal  powers, 
frequencies,  arrival  angles,  and  bandwidths.  The  results  will  be  of 
interest  to  both  the  theoretician  and  the  system  designer. 

We  use  P(e),  the  probability  of  a  symbol  error  at  the  receiver 
output,  as  the  measure  of  system  performance.  For  each  combination  of 
desired  signal  and  interference,  we  calculate  the  signals  at  the  array 
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output.  We  then  calculate  P(e)  of  the  detector  when  the  array  output 
signals  appear  at  its  input.  In  each  case  we  assume  that  the  system 
uses  the  ideal  detector  for  the  desired  signal  in  the  presence  of  white 
Gassian  noise. 

Although  a  high  SINR  at  the  array  output  is  generally  desirable,  we 
show  that  P(e)  at  the  detector  output  cannot  generally  be  determined 
from  the  SINR  alone.  For  example,  the  performance  of  most  digital 
detectors  with  Gaussian  noise  and  continuous  wave  (CW)  interference 
depends  on  several  factors  including  desired  signal  power,  interference 
power,  noise  spectral  density,  and  noise  bandwidth.  If  we  hold  all  but 
one  of  these  parameters  constant  and  let  the  remaining  parameter  vary, 
then  we  find  that  P(e)  does  vary  monotonical ly  with  SINR.  However,  if 
we  allow  more  than  one  parameter  to  vary  at  a  time,  we  find  it 
impossible  to  predict  P(e)  from  the  SINR  alone.  In  fact,  there  are  many 
cases  where  P(e)  becomes  worse  as  the  SINR  increases.  Therefore,  we  use 
P(e)  and  not  SINR  as  the  measure  of  performance. 

We  find  that  the  performance  of  a  communication  system  that 
includes  an  LMS  array  is  strongly  dependent  upon  the  bandwidth  of  the 
signals  at  the  array  input.  Generally,  we  find  the  system  performance 
improves  as  the  array  bandwidth  decreases  (as  long  as  this  bandwidth  is 
large  enough  to  pass  the  desired  signal).  However,  many  systems  require 
an  array  bandwidth  that  is  several  times  larger  than  the  data  bandwidth 
of  the  transmitted  bit  stream.  There  are  many  reasons  why  this  might  be 
the  case.  First,  the  spectrum  of  a  signal  modulated  by  any  of  the 
classical  digital  modulation  techniques  discussed  above  is  not  limited 
to  the  data  bandwidth.  Such  a  signal  will  be  distorted  if  it  is 
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filtered  to  occupy  a  narrow  bandwidth.  Second,  it  may  be  difficult  to 
implement  the  filters  required  at  the  array  input  to  limit  the  array 
bandwidth.  Third,  the  system  may  be  required  to  operate  over  a  wide 
frequency  range  to  accommodate  several  signal  channels  or  unknown 
doppler  shifts.  Finally,  the  spectrum  of  the  desired  signal  may  be 
intentionally  spread  for  data  security  or  to  permit  the  generation  of  a 
reference  signal  [3,4,13].  For  these  reasons  we  examine  the  system 
performance  for  several  array  bandwidths. 

We  first  examine  the  performance  of  wideband  BPSK,  QPSK,  OPSK,  and 
FSK  systems  with  CW  interference.  Ry  wideband  we  mean  that  the  array 
bandwidth  in  Hz  is  at  least  ten  times  as  large  as  the  transmitted  bit 
rate  in  bits/sec.  We  find  that  the  qualitative  performance  of  the 
various  systems  are  commensurate  with  the  well-known  performance 
characteristics  of  their  respective  modulation  techniques.  Therefore, 
for  the  remainder  of  the  study,  we  use  BPSK  modulation  with  the 
understanding  that  similar  results  would  be  obtained  using  other 
modulation  methods. 

We  examine  several  array  bandwidths  using  the  wideband  models.  In 
each  case  we  find  that  the  system  performance  improves  as  we  reduce  the 
array  bandwidth.  In  order  to  determine  the  optimum  array  bandwidth  we 
develop  a  narrowband  BPSK  modulation  method.  This  method  allows  the 
desired  signal  to  occupy  a  channel  with  a  bandwidth  equal  to  the  desired 
signal  data  rate.  We  examine  the  system  performance  for  a  wide  range  of 
array  input  bandwidths.  Under  most  interference  scenarios,  we  find  the 
best  performance  when  the  array  bandwidth  is  as  small  as  possible  (i.e.. 
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when  the  array  bandwidth  is  equal  to  the  data  bandwidth).  However,  for 
weak  interference,  this  general  rule  does  not  always  hold.  We  examine 
the  reasons  for  this  anomalous  behavior. 

In  addition  to  CW  interference,  we  also  examine  the  effects  of 
bandlimited  noise  interference  on  the  RPSK  system  performance.  We  find 
that  the  LMS  array  with  a  single  complex  weight  behind  each  element 
often  has  diffiulty  nulling  strong  broadband  interference.  This  problem 
is  overcome  if  we  use  tapped  delay  line  filters  to  implement  the  weights 
[15].  We  compare  the  performance  of  arrays  using  both  of  these  weight 
implementations  for  various  signals  powers,  bandwidths,  and  arrival 
angles. 

Chapter  II  presents  a  description  of  the  LMS  array  and  the 
calculations  required  to  predict  the  signal  power  levels  at  the  array 
output.  Chapters  III  through  VI  describe  the  performance  of  the  RPSK, 
QPSK,  OPSK,  and  FSK  systems  with  CW  interference.  Chapter  VII  examines 
the  performance  of  a  narrowband  RPSK  system  with  CW  interference  for 
various  input  bandwidths.  Chapter  VIII  describes  the  performance  of  a 
RPSK  system  with  broadband  interference  for  both  the  LMS  and  the  tapped 
delay  line  LMS  arrays.  Finally,  Chapter  IX  contains  the  conclusions. 
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CHAPTER  II 


THE  LMS  ADAPTIVE  ARRAY 


A.  INTRODUCTION 

In  this  chapter  we  will  describe  the  basic  operating  principles  of 
the  LMS  adaptive  array.  We  present  the  mathematical  models  used  in 
later  sections  to  predict  array  performance. 

The  discussion  of  array  operation  presented  in  this  chapter  is 
abbreviated  and  is  presented  mainly  to  define  the  notation  used  in  the 
following  chapters.  More  detailed  descriptions  of  the  LMS  array  and 
its  operation  are  given  in  [2]  and  [15].  The  notation  used  in  this 
section  is  generally  consistent  with  that  used  in  these  references. 

B.  THE  COMPLEX  LMS  ARRAY 

Figure  2.1  shows  a  block  diagram  of  an  N-element  adaptive  array. 
The  signal  from  each  antenna  element  is  fed  through  an  ideal  bandpass 
filter  (BPF),  which  limits  the  thermal  noise  power  and  rejects  out-of- 
band  interference  at  the  array  input.  We  denote  the  analytic  signal  at 
the  input  to  the  ith  array  input  by  x^t).  For  the  N-element  array  we 
define  the  signal  vector  X  which  is  composed  of  the  array  input 
signals. 
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Figure  2.1.  The  complex  IMS  array. 
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Xj(t) 

x2(t) 

_  *N(t)  __ 


(2.1) 


Each  element  signal  is  multiplied  by  a  complex  weight  and  the  weighted 
signals  are  summed  to  produce  the  array  output  s(t).  We  define  the 
complex  weight  vector. 


W  = 


w 


2 


(2.2) 


The  array  output  is  given  by, 

s(t)  =  WTX  .  (2.3) 

The  array  output  is  subtracted  from  a  reference  signal,  r(t),  to 
produce  an  error  signal  e(t), 

e(t)  =  r(t)  -  s(t)  .  (2.4) 


The  IMS  array  adjusts  the  array  weights  to  minimize  the  mean  square 
value  of  this  error  signal  [2].  A  suitable  reference  signal  for  the 
array  can  often  be  extracted  from  the  signals  at  the  array  output.  The 
circuit  that  derives  this  reference  signal  typically  uses  known 
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properties  of  the  desired  signal  (such  as  spread-spectrum  sequence 
characteristics)  to  partially  separate  the  desired  and  undesired 
components  of  the  array  output  [3,4,5].  The  signal  at  the  circuit 
output  should  be  well  correlated  with  the  desired  signal  and  poorly 
correlated  with  the  interference.  The  references  cited  in  Chapter  I 
discuss  reference  generation  techniques  for  several  types  of  desired 
signal  modulation.  In  this  report  we  do  not  address  the  problem  of 
reference  signal  generation.  We  simply  assume  that  the  reference  signal 
is  a  perfect  replica  of  the  desired  signal  at  the  output  of  BFP1. 

In  steady  state  the  LMS  weight  vector  is  given  by  [1], 

W  =  *_1S  ,  (2.5) 

where  >5  is  the  Covariance  Matrix, 

5  =  E[X*xT]  ,  (2.6) 

and  S  is  the  reference  correlation  vector,  which  is  given  by, 

S  =  E[X*r(t)]  .  (2.7) 

where  E[®]  denotes  the  expected  value.  It  can  be  shown  [15]  that  the 
steady  state  weights  given  by  (2.5)  produce  the  maximum  attainable 
Signal -to-Interference-pl us-Noi se-Ratio  (SINR)  at  the  array  output. 

It  is  often  convenient  to  separate  the  desired,  interfering,  and 
thermal  noise  components  of  the  incoming  signal  vector  as, 

x  *  xo  *  xi  *  *n  <2-s> 
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If  we  assume  that  the  desired,  interfering,  and  thermal  noise  signals 
are  uncorrelated  (and  zero  mean),  then  we  can  separate  the  three  related 
components  of  the  covariance  matrix 

*  =  E[xV] 

-  E[Xd*XdT]  *  ECx/x,1]  ♦  ECX„*XnT7  . 

■  ‘d  +  *1  *  4n 

where 

5d  =  E[Xd*XdT]  1  (2An) 

»  E[X1*XiT]  ,  (2.11) 

and 

%  =  E[Xn*XnT]  .  (ZAP.) 

Thus  we  have  divided  the  covariance  matrix  into  a  sum  of  three  separate 
matrices,  each  involving  only  one  signal.  This  step  simplifies 
subsequent  derivations. 

We  assume  in  the  derivations  below  that  we  have  a  three 
element  linear  array  with  one  half  wavelength  spacing  between  elements. 
Figure  2.1  shows  the  array  geometry.  Of  course,  we  can  extend  our 
results  to  cases  where  the  array  has  more  than  three  elements;  however, 
the  three  element  array  calculations  suitably  illustrate  the  system 
operating  principles  and  performance. 

We  assume  in  the  derivations  below  that  we  have  a  3-element  linear 
array  with  one-half  wavelength  element  spacing  as  shown  in  Figure  2.1. 
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A  3-element  array  has  two  degrees  of  freedom.  The  array  can  therefore 
steer  a  null  on  the  interference  and  simultaneously  maximize  the  desired 
signal  response.* 

We  assume  that  the  desired  signal  is  incident  from  an  angle  of  0d 
(measured  from  the  broadside  direction)  and  a  single  interfering  signal 
is  incident  from  an  angle  of  0^. 

The  interelement  time  delay  between  the  signals  on 
adjacent  elements  are  given  by 

Td  =  (L/c)sined  (2.13) 

and , 

T.  =  (L/c)sinei  (2.14) 

for  the  desired  signal  and  interference  respectively,  where  c  is  the 
propagation  velocity  and  l  is  the  element  spacing.  The  interelement 
phase  shifts  for  the  desired  and  interfering  signals  are  given  by, 

♦d  =  1TSln(0d)  (2.16) 

and , 

=  1I((JJ-j/wd)sin(9i )  (2.16) 

respectively,  where  are  the  desired  signal  and  interference 

frequencies.  Slight  variations  in  uij  about  u>d  have  little  effect  on 

♦There  are  limits  to  the  resolution  capabilities  of  the  3-element  array. 
We  shall  find  in  later  chapters  that  the  array  cannot  keep  the  desired 
signal  at  a  pattern  maximum  when  a  very  strong  interfering  signal  is 
spatially  close  to  the  desired  signal. 
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most  performance  calculations.  From  Equation  (2.16)  we  see  that  a  small 
change  in  oj,  produces  a  corresponding  change  in  <j>j.  Similarly  a  small 
change  in  0,  has  the  same  effect.  Therefore,  the  change  in  the  array 
performance  caused  by  variations  in  u>j  are  identical  to  those  caused  by 
a  corresponding  variation  in  9^.  If  these  variations  are  small,  then 
the  array  performance  is  negligibly  affected.  In  the  following  chapters 
we  usually  assume  that  <t>i=7rsin0-j .  This  assumption  will  allow  us  to 
calculate  the  system  performance  without  specifying  the  operating 
frequency. 

We  denote  the  analytic  signal  representation  of  the  desired,  and 
interfering  signals  as  received  by  element  1  by  3(t)  and  T(t) 
respectively.  We  define  the  complex  autocorrelation  function  (ACF)  for 
the  desired  and  interfering  signals  as, 

Rd(0  =  E[d*(t)d(t+x)]  (2.17) 

and , 

R?(t)  =  E[T*(t)T(t+x)]  .  (2.18) 

We  assume  that  white  Gaussian  noise  is  present  at  the  input  of  each  of 
the  element  BPF's.  We  assume  that  the  noise  signals  at  the  various  8PF 
outputs  are  mutually  statistically  independent.  We  denote  these  signals 
by  n^(t),  ^(t)  and  ^(t)  with, 

E[Ki(t)n.(t)]  =  o2sij  ,  (2. IP) 

where  5^  is  the  Kronecker  delta. 
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We  now  can  write  the  desired  signal,  interference  signal,  and  noise 


covariance  matrices  for  the  three-element  array. 


-  Rd(0) 

Rdt-Td) 

Rd(-2Td)  - 

id  = 

Rd(Td) 

«d(0) 

Rd(-Td) 

_  Rd(2Td) 

Rd(Td) 

Rd(0) 

-  rHo) 

RT(-fi) 

R?(-2Ti)  - 

*i  - 

R?(Ti) 

«T(0) 

R?(-Td) 

R?(2Ti) 

R?(0) 

(2.20) 


(2.21) 


and. 


where  I  is  the  identity  matrix. 

The  covariance  matrix  5  can  now  be  determined  from  (2.9),  and 
(2.20)  through  (2.22). 

We  assume  the  reference  signal  for  the  array  is  of  the  form, 


r(t)  =  R  d(t) 


(2.23) 


We  note  that  the  SINR  at  the  array  output  is  independent  of  the 
magnitude  of  R  (for  non-zero  R).  An  increase  or  decrease  in  R  produces 
a  proportionate  increase  or  decrease  in  the  level  of  each  signal  at  the 
array  output.  We  are  therefore  free  to  arbitrarily  set  R  without 
affecting  the  relative  values  of  the  signals  at  the  array  output. 
Typically  a  reference  generation  system  generates  a  reference  signal 
which  is  highly  correlated  with  the  desired  signal  and  has  fixed 
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amplitude.  We  can  model  such  a  fixed-amplitude  signal  by  letting 


R  =  1// ( Pd ) i n  where  (pd)in  is  the  power  of  the  desired  signal. 

With  the  reference  signal  of  the  form  given  by  (2.23)  the  reference 
correlation  vector  is  given  by, 

S  =  E[X*r(t)] 

=  R  E[X*d(t)]  .  (2. P4) 

Since  d(t),  T(t),  and  the  noise  signals  are  uncorrelated  we  can  expand 

(2.24)  into, 

S  =  R  E[(Xd*  +  Xi*  +  Xn*)d(t)] 

«  R  E[Xd*d(t)]  .  (2.25) 

Thus,  for  our  three  element  array 

-Rd(o)  - 

S  =  R  Rd(Td)  (2.25) 

_  Rd(2Td)  _  . 

The  complex  weight  vector  is  then  determined  from  (2.5),  (2.9),  and 

(2.25) . 

We  can  now  determine  the  desired,  interfering  and  noise  signals  at 
the  array  output.  From  (2.3)  and  (2.8)  we  write  the  array  output  signal 
s(t)  as, 

5(t)  =  sd(t)  +  s^(t)  +  s„(t)  (2.27) 
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where. 


3 

Sj(t)  =  I  w  3[t-(j-l)T ,]  (2.28) 

j=l  J 


3 

S?(t)  -  l  w  T[t-(j-l)Ti]  (2.29) 

J=1  J 

and , 


Sq(t)  =  wjni(t)  +  W2h2(t)  +  w3n3(t)  .  (2.30) 


The  power  of  the  desired  signal  at  the  array  output  is  given  by. 


Pd  =  l/2Rd(0) 


»  l/2E{[sd(t)][sd(t)]*}  .  (2.31) 


From  (2.28)  and  (2.31),  we  have 


3  3 

Pd  =  1/2  l  l  w  .wk*Rd[(j-k)T  )  (2.32) 

j-1  k=l 


Similarly  we  have. 


3  3 

P.  =1/2  l  l  wiw.*R~[(j-k)Ti1  (2.33) 

j-1  k=l 


and , 


P 


n 


3 

1/2  l 
j=l 


,2  2 

"j1  5 


(2.34) 


15 


The  SINR  at  the  array  output  is  given  by 


SI  NR  = 


P 


d 


pi  +  Pn 


(2.35) 


In  many  of  the  previous  adaptive  array  studies,  the  SINR  at  the 
array  output  was  used  as  a  figure  of  merit.  In  the  following  sections 
we  examine  several  different  types  of  digital  communication  systems  and 
determine  their  performance  using  the  probability  of  error  as  a  figure 
of  merit.  We  show  that,  although  SINR  is  often  useful,  it  does  not 
always  give  sufficient  information  to  accurately  predict  the  system 
performance. 
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CHAPTER  III 


PERFORMANCE  OF  A  BPSK  SYSTEM  WITH  CM  JAMMING 

A.  INTRODUCTION 

In  this  chapter  we  examine  the  performance  of  a  BPSK  communication 
system  that  includes  an  LMS  array  to  suppress  CW  interference.  Figure 
3.1(a)  shows  a  typical  BPSK  signal  waveform.  A  "1“  symbol  is 
represented  by  a  sinusoidal  waveform  of  duration  T  seconds.  A  "0" 
symbol  is  represented  by  a  similar  sinusoid  waveform  which  is  180°  out 
of  phase  with  the  "1"  waveform.  The  receiver  estimates  the  phase  of  the 
received  signal  at  the  end  of  each  symbol  interval  and  decides  whether  a 
"1"  or  "0"  was  sent  during  the  interval. 

The  model  we  develop  in  this  chapter  assumes  that  the  bandwidths  of 
the  BPFs  at  the  element  inputs  are  several  times  larger  than  the  data 
bandwidth  of  the  transmitted  bit  stream.  There  are  several  reasons  why 
this  assumption  accurately  models  many  real  systems.  First,  the 
spectrum  of  a  BPSK  signal  is  not  limited  to  the  data  bandwidth.  Such  a 
signal  will  be  distorted  if  it  is  filtered  to  occupy  a  narrow  bandwidth. 
Second,  it  may  be  difficult  to  implement  the  filters  required  at  the 
array  input  to  limit  the  array  bandwidth.  Third,  the  system  may  be 
required  to  operate  over  a  wide  frequency  range  to  accommodate  several 
signal  channels  or  unknown  doppler  shifts.  Finally,  the  spectrum  of  the 
desired  signal  may  be  intentionally  spread  for  data  security  [161  or  to 
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permit  the  generation  of  a  reference  signal.  In  the  last  case  the 
spectrum  may  be  spread  by  either  di rect -sequence  spreading  [3,4]  or 
frequency  hopping  [14].  We  defer  consideration  of  more  narrowband 
systems  (those  with  input  bandiwdths  that  are  approximately  equal  to  the 
data  bandwidth)  until  Chapter  VII. 

In  Section  8  we  examine  the  performance  of  the  LMS  array  with  the 
desired  and  interfering  signals  as  described  above.  In  Section  C  we 
calculate  P(e)  for  the  ideal  BPSK  detector  when  the  desired  signal  is 
corrupted  by  both  CW  interference  and  Gaussian  noise.  In  Section  n  we 
determine  the  performance  of  the  combined  LMS  array  and  BPSK  detector. 

We  determine  how  this  performance  depends  on  the  received  signal  power 
levels,  frequencies,  and  the  bit  duration. 


B.  ARRAY  PERFORMANCE 

During  the  nth  symbol  interval  the  BPSK  desired  signal  at  the 
output  of  BPF1  is  given  by, 

d(t)  =  Adexp{j[u)dt  +  <J>n(t)  +  <i>d]}  (n-l)T<t<nT  (3.1) 

where 

Ad  =  the  desired  signal  amplitude  at  the  array  input, 

=  the  desired  signal  carrier  frequency, 

ipd  =  a  random  phase  angle  uniformly  distributed  on  [0,2*] 

(We  denote  this  by  saying  Ud  is  U[0,2*].) 

4>n(t)  =  the  data  modulation  during  the  n^  symbol  interval  where, 

,  0  with  prob.  1/2  (during  "0"  symbols) 

»nu)  -  i 

*  with  prob.  1/2  (during  "1"  symbols). 
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and  T  is  the  symbol  duration. 

We  assume  that  the  symbol  stream  is  stationary.  The  complex 
autocorrelation  function  for  the  desired  signal  is  given  by, 

Rj(x)  =  E[d*(T)d(t+x)] 

’  Ad2[l  -  |t|/T]  exp[ju)dT]  | t |  <  T 

(3.?) 

,0  1  r  |  >  0  . 

The  Power  Spectral  Density  (PSD)  of  the  desired  signal  is  shown  in 
Figure  3.1(b). 

i(t),  the  CW  interfering  signal,  is  given  by, 

?(t)  =  A.  exp[j<jJit+’4>i 3  (3.3) 

where, 

=  the  interference  amplitude, 

tp .  *  the  interference  phase.  is  U[0,2tt]  and 
1  independent  of  <j>d. 

and , 

u)_.  =  the  interference  frequency. 

The  complex  autocorrelation  function  of  this  signal  is  given  by, 

R?(t)  =  E[T*(t)1(t+r)] 

=  Ai2exp[ju)^T3  . 


(3.4) 


We  assume  an  independent  thermal  noise  component  with  power  a2  is 
present  in  each  element  signal.  We  therefore  have 


E[ni*(t)nJ.(t)]  =  o25ij.  . 


(3.5) 


We  now  calculate  the  array  weights  and  the  power  levels  of  the 

*N» 

various  signals  at  the  array  output.  Since  d(t),  i(t)  and  the  thermal 
noise  are  uncorrelated  we  can  determine  the  covariance  matrix  by  summing 
its  three  components. 


[Id|,  j$d 


'd  d 


A/  I  (l-T-)® 


, .  ilii  -j*d 

(1-  j  )e 


(l_]2Tdl)eJ2,d 


(3.6) 


7  Ji’i 
9i  =  Ai2  e 


(3.7) 


J4>  i 

e  e 


v*2' 


(3.S) 


where  I  is  the  identity  matrix  and  it  is  assumed  that  1 2Tci | <T . 


We  make  the  simplifying  assumption  that  the  symbol  period  is  very 
long  compared  with  the  interelement  time  delay.  This  assumption  is 
reasonable  since  the  elements  are  one-half  wavelength  apart.  Thus  only 
a  fraction  of  an  RF  carrier  cycle  occurs  during  a  time  interval  of 
duration  T^.  Many  such  cycles  occur,  however,  during  each  bit  interval 
of  length  T.  With  this  assumption  we  can  calculate  the  array 
performance  without  specifying  the  operating  frequency  or  symbol  rate. 
The  desired  signal  covariance  matrix  becomes 


-j*d 

e 

-j  2<}><j 
e 

j<t>d 

1 

-Hd 

e 

J  2<i>d 

j<t>d 

e 

1 

(3.0) 


The  total  covariance  matrix  then  becomes. 


+  ■#_ 
din 


A  .2  +A . 2  +o2 
d  i 

-  y  y  -J^i 

A  2e  +A.2e 

d  i 

A  2e*j2*d+A.2e"J?>i 

d  i 

.  y  j^d  .  y  •i'H 

A  2e  +A.2e 

A  ,2+A.2+o2 

A  2e’j,t>()  +  A.2_jH 

d  i 

d  i 

d  i 

VeJ2*«.Vej2*’ 

J^d  .  y  3 i 

A  2  e  +A . 2  e 

d  i 

A,2  +  A.2  +  o2 
d  i  _ 

(3.10) 
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The  reference  correlation  vector  (2.2fi)  is  given  by. 


“  1 


exp[j<j.d] 

exp[j24>d] 


(3.11) 


We  can  now  calculate  the  weight  vector  from  (2.5)  and  Pd,  p^  ,  and  Pn 
from  (2.32),  (2.33)  and  (2.34). 

In  the  next  section  we  derive  the  performance  of  the  ideal  BPSK 
detector  when  the  desired  signal  is  corrupted  by  thermal  noise  and  CW 
interference.  In  order  to  determine  the  detector  performance  we  must 
know  the  composite  signal  at  the  detector  input.  To  determine  this 
signal  we  switch  to  the  real  representation  of  the  signals  at  the  array 
output. 

At  the  array  output  the  desired  signal  is  a  BPSK  signal  with  power 
Pd .  This  (real)  signal  during  the  nt*1  bit  interval  is  given  by. 


sd(t)  =  /Wd  cosOdt+$n(t)+Ydl» 


(3.12) 


where  Pd,  <*>d,  and  $n(t)  are  as  defined  above  and  yd  is  the  desired 
signal  phase  angle  at  the  array  output.  Since  the  covariance  matrix 
and  the  reference  correlation  vector  are  both  independent  of  4>d  (the 
phase  angle  of  the  desired  signal  at  the  array  input),  the  steady  state 
weights  are  also  independent  of  !>d.  Therefore,  yd  is  given  by  (yd)0  + 

■J>d  where  (yd)0  is  some  fixed  phase  angle.  Thus,  since  i|>d  is  U[0,2^t], 
the  (modulo-2it)  value  of  yd  is  also  !J[0,2ir],  This  is  a  property  of  the 
modulo-M  sum  of  two  variables,  one  of  which  is  1J[0,M],  This  property  is 
proven  in  [17]  by  Scire. 
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We  now  determine  the  amplitude  and  phase  of  the  CW  interfering 
signal  at  the  array  output.  This  signal,  which  has  power  P^ ,  is  given 

by. 


Si(t)  =  /2PiCOs[>i  (t)  +  Yi  U  .  (3.13) 

where  y-j  is  the  phase  angle  of  the  interference  at  the  array  output.  We 
can  show  that  y^  is  U[0,2?r]  using  arguments  analogous  to  those  used  above 
to  show  that  y<j  is  U[Q,2ir]. 

Finally,  the  noise  signal  at  the  array  output  is  a  zero-mean 
Gaussian  random  process,  sn(t),  with  variance  Pn.  This  noise  is 
bandlimited  Gaussian  noise  with  a  flat  PSD  in  the  passband  of  the  ideal 
BPFs  at  the  element  inputs.  The  value  of  the  two-sided  PSD  in  this 
passband,  which  we  denote  by  n/2  is  given  by, 

N0  3 

n/2  =  l  |wm|2  (3.14) 

m=l 

where  N0/2  is  the  two-sided  thermal  noise  PSD  at  the  element  inputs.  The 
PSD  of  the  noise  signal  at  the  array  output  is  shown  in  Figure  3.2(a). 

C.  BPSK  DETECTOR  PERFORMANCE  WITH  CW  AND  NOISE  INTERFERENCE 

In  this  section  we  derive  the  performance  of  an  ideal  BPSK  detector 
when  the  BPSK  desired  signal  is  corrupted  by  additive  white  Gaussian 
noise  (AWGN)  and  CW  interference.  Rosenbaum  [IB]  and  Prabhu  [IB]  have 
examined  similar  problems;  however,  these  early  works  were  typically 
rather  vague  in  their  description  of  the  "ideal"  receiver.  The  receiver 
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model  in  these  papers  assumes  that,  in  addition  to  a  desired  and 
interfering  signals,  there  is  Gaussian  noise  with  power  a2  at  the 
detector  input.  However,  the  importance  of  the  noise  PSD  and  the  signal 
energy-per-bi t  are  generally  not  discussed  in  detail.  The  dependence  of 
performance  on  the  frequency  of  the  interfering  signal  is  also  not 
di scussed. 

To  thoroughly  understand  the  effects  of  the  various  signal 
energies,  the  noise  PSD,  and  the  interference  frequency  on  the  detector 
performance  we  shall  carefully  derive  P(e)  of  the  ideal  RPSK  detector. 
The  RPSK  detector  shown  in  Figure  3.3  is  ideal  if  the  only  undesired 
signal  is  additive  white  Gaussian  noise  (AWGN)  [20],  The  incoming 
signal,  s(t),  is  multiplied  by  cos  (u^t+Yd)  •  The  output  of  the 
multiplier  is  integrated  over  the  bit  interval  and  a  bit  decision  is 
made  based  upon  the  sign  of  the  integrator  output  at  the  end  of  the  bit 
interval.  A  positive  output  indicates  a  "0"  bit  and  a  negative  output 
indicates  a  "1"  bit. 

Note  that  this  ideal  detector  requires  a  coherent  reference  signal 
to  operate.  A  method  for  deriving  such  a  signal  from  the  received 
signal  is  discussed  in  [3].  We  assume  that  a  suitable  reference  signal 
is  available  and  that  this  signal  is  free  of  amplitude  and  phase  noise. 

Summarizing  the  results  of  the  previous  sections,  the  composite 
signal  at  the  detector  input  is 

s(t)  =  sd(t)  +  si (t)  +  sn(t) 


2S 


(3.1b) 


(a)  Noise  spectrum  at  detector  input. 


PSD  (W/Hz) 


(b)  Noise  spectrum  at  integrator  input. 


(c)  I H (u> ) ! 2  for  ideal  integrator. 

Figure  3.2.  Noise  spectrum  and  transfer  function  for  ideal  integrator. 
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where  sd(t)  and  s^(t)  are  given  by  (3.12)  and  (3.13),  and  sn(t)  is  a 
narrowband  Gaussian  random  process  which  has  a  flat  (two-sided)  PSD  with 
magnitude  n/2  in  the  passband. 

The  signal  at  the  input  to  the  integrator  is  given  by 

S(t)C0S(a)dt+Yd)  =  [sd(t)+s1  (t)+Sn(t)]  COS(udt+Yd)  (3.1b) 

We  can  substitute  sd(t)  and  s,(t)  into  this  expression  and 
integrate  between  (n-l)T  and  nT  to  determine  the  signal  at  the 
integrator  output  at  the  end  of  the  bit  interval 

S(n)  =  Xd(n)+Xi(n)+Xn(n)  (3.17) 


where. 


Xd(n) 


nT  _ 

/  /2Pdcos[u>dt  +  4>n (t )  +  Ydl  cos[u)dt  +  Yd]  dt 

(n-l)T 


(3.18) 


nT  _ 

X.(n)  =  /  v'2Pi  cos[a)it  +  y i ]  cos[a!dt  +  Yd]  dt 
(n-l)T 


and. 


nT 

xn(n)  =  /  sn(t)cos[u)dt+Yd]  dt 
(n-l)T 


(3.19) 


(3.20) 
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Figure  3.3.  The  ideal  BPSK  detector. 


We  now  evaluate  Xd(n),  X-j(n)  and  Xn(n)*  and  then  combine  the  results  to 
determine  X(n).  Application  of  a  standard  trigonometric  identify  to  the 
integrand  in  Equation  (3.18)  yields, 

_  nT 

xd(n)  =  ^Pd/2  /  cosf>dt+4>n(t)  +  2yd]  +  cos[*  (t)1  dt  .  (3.21) 

(n-l)T 


*  Note  that  Xd(n),  X-j(n)  and  Xn(n)  are  not  the  signal  vectors  Xd,  X,  and 

Xn  defined  in  Chapter  2.  Although  the  similar  notation  might  be 
somewhat  confusing  at  first,  the  intended  variable  should  be  clear 
from  the  context. 
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We  ignore  the  double  frequency  term  since  its  contribution  to  the 
value  of  the  integral  is  negligible  due  to  our  assumption  that  many 
carrier  cycles  occur  during  each  bit  interval.  Therefore,  we  have, 

nT 

xd(n)  =  /pd/2  I  cos[*n(t)]  dt 

(n-l)T 

=  /T^Tzjl  cos[>n(t)]  ,  (3.22) 

since  4>  ( t )  is  constant  during  the  bit  interval. 

We  now  calculate  X^(n).  We  first  expand  the  integrand  in  Equation 
(3.19), 

nT 

xi  (n)  =  /P./2  f  {cos[(u)d+u).)t  +  Yd  +  y.] 

(n-l)T 

+  COS[(u)d-uK)t  +  Yd  -  Ti  1}  dt  .  (3.23) 

We  let 


Yrel  Yd  "  Ti  ’ 


(3.24) 


and , 


Ao) 


"  “i  * 


(3.2b) 
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Substitution  of  these  new  variables  into  (3.23)  yields, 

_  nT 

X  (n)  =  / P- / 2  /  cos ( Aujt  +Yre1)  dt.  (3.26) 

1  (n-l)T 

after  we  drop  the  second  harmonic  terms.  We  evaluate  the  integral  to 
give  us, 

Xi(n)  -  [sin(nAu,T  +  yre1 )  -  Sin(n-l)Aa>T  +  Yrel  ]  (3.27) 

Simplification  of  this  expression  with  standard  trigonometric  identities 
yields, 

___ „  AujT  2n-l 

Xi  (n)  -  sine  (— )  cos  (yrel  +  —  AcuT),  (3.23) 

where  the  sine  function  is  given  by, 

sinc(x)  =  [sin(x)]/x,  (3. 2d) 

with 

sinc(O)  =  1.  (3.30) 

From  (3.28)  we  see  that,  for  a  given  interference  power  level,  the 
amplitude  of  the  interference  at  the  integrator  output  is  directly 
proportional  to  sinc(Au>T/2).  At  the  points  where  the  sine  function  is 
zero  (i.e.  for  AooT/2  =  nir)  the  interference  does  not  disrupt  the 
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detector.  At  these  values  of  AwT,  the  baseband  interference  product  at 
the  multiplier  output  is  a  sinusoid  with  an  integral  number  of  cycles 
occurring  during  the  integration  period.  Therefore  the  contribution  of 
the  interference  to  the  integrator  output  is  zero.  The  third  integral, 

X  (n),- is  the  integral  of  the  product  sn (Ocosfu^t+y^) .  X^(n)  is  most 

easily  evaluated  in  the  frequency  domain.  If  the  input  noise  signal  has 
the  PSD  shown  in  Figure  3.2(a)  then  the  signal  at  the  integrator  input 
has  the  PSD  N(w)  shown  in  Figure  3.2(b). 

In  the  frequency  domain  the  integrator  is  equivalent  to  a  linear 
filter  with  a  transfer  function  of 


H(u) 


l-exp(jaiT) 


(3.31) 


At  the  end  of  the  integration  period  the  noise  produces  an  additive 
Gaussian  random  variable  at  the  integrator  output  This  zero  mean  random 
variable  has  a  variance  N',  which  is  given  by. 


N'  =  (^)  J  N(oi)|H(a.)|2  dcu 


1  oo  U)X 

=  (“ ■)  T2  /  N(w)sinc2(^~)  dio.  (3.32) 


We  can  make  some  simplifying  assumptions  which  make  evaluation  of 

(3.32)  more  strai ghtforward  with  only  negligible  effects  on  the  results. 

OJT 

In  Figure  3.2(c)  we  show  sinc2(^~).  The  integrand  in  (3.32)  is  the 
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product  of  the  two  curves  shown  in  Figure  3.2(b)  and  3.2(c).  From 

ojT 

these  figures  it  is  apparent  that  the  contribution  of  N(uo)si nc2  (~)  to 
the  value  of  the  integral  quickly  diminishes  as  u>  becomes  large  (due  to 
the  fast  roll-off  of  the  sine2  term).  Since  we  have  assumed  that  the 
input  BPFs  are  much  larger  than  the  data  bandwidth  of  the  desired 

signal,  we  can  make  the  approximation  that  the  noise  at  the  integrator 

n 

input  is  white,  with  PSD  —  at  all  frequencies.  With  this  approximation 
we  can  easily  evaluate  (3.32) 


1  00  T 

N'  =  (^*)T2(^)  /  sinc(^r)  d«. 
-00 

If  we  now  make  the  substitution, 
wT 


(3.33) 


(3.34) 


then  (3.33)  becomes, 

N'  =  (7)(J)T  /  si nc2  (x)  dx.  (3.35) 

-00 

The  integral  in  (3.35)  is  found  to  be  equal  to  tt  using  Parseval's 
theorem  [21].  Thus  N',  the  variance  of  the  noise  random  variable  at  the 
integrator  output,  is  given  by 


N'  =  Tn/4.  (3.36) 

The  error  in  N'  due  to  our  assumption  that  the  noise  at  the  detector 
input  is  white  is  less  than  2%  if  the  actual  noise  bandwidth  at  this 
input  is  greater  than  approximately  10/T  Hz.  The  effects  of  noise 
bandwidth  are  discussed  further  in  the  next  section. 
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We  now  combine  (3.22),  (3.28),  and  (3.36)  to  form  the  expression 
for  the  composite  random  variable  at  the  integrator  output  at  the 
sampling  instant.  This  random  variable,  X(n),  is  given  by, 

X(n)  =  /TP72)Tcos(>n(t)l 

?n  1 

+  /(Pi/2)Tsinc(Aa1T/2)cos(Yrel  +  — AwT  +  Xp(n) ,  (3.37) 

where  1>n(t)  is  either  0  or  t  depending  on  the  data  bit  and  X  (n)  is 
normally  distributed  with  zero  mean  and  variance  nT/4. 

We  now  can  determine  the  probability  of  error  for  the  ideal  BPSK 
receiver.  If  both  bits  are  equally  likely  the  conditional  probability 
of  error  given  y^  is  given  by, 

P(e|yrel)  =  0.5P(e|yre1,  “0"  sent)  +  0.5P(e|Yrel ,  "1"  sent)  . 

=  O.^P{X(n)  <  0|Yre1,*n(t)«0}  +  O.SPt X(n)>0(yre1  ,4>n(t)=^} 

P(e|Yrel)  =  0.5P(/lP^72lT  + /TV2)Tsinc(Aa,T/2)cos(yrel+ 

+  Xn(n)  <  0} 

Pp  1 

+  0. 5P{  V(Pd/2)T  +  /(P1/2)Tsinc(AuJT/2)cos(yrel  +  ~ — WT) 

+  Xn(n)  >  0}. 

p(e|Yrel)  -  0.5P(Xn(n)  <  - 

, _  2n-l 

-  »  (Pi /2)Tsinc(AwT/2)cos(yre1  +  AuT)} 

+  °- 5Pt  xn(n)  >  ST vF)t 

_  2n-l 

-  / (P.  /2)Tsinc(Aa)T/2)cos(yrel  +  —  Au>T)} 
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P(e|Yrei)  =  °«5erfc 


,A(oT.  .  2n-l 

nc(— )cos(yrei+  ~Y~ AuT) 


+  0.5  erfc 


,  /  2PdT  J^V  .  .  2n-l  . 

V  T~  +  V  ~ n  sinc(Au)T)cos(Yre^  +^~ Au)T) 


where  erfc  is  the  complementary  error  function  defined  by, 
erfc(x)  =V(27)  /  exp(-z2/2)  dz. 


(3.38) 


(3.39) 


We  now  define  the  desired  signal  and  interference  energy-per-bit- 
interval  by. 


Ed  *  PdT* 


and 


Ei  *  V- 

respectively.  Equation  (3.38)  can  then  be  rewritten  as. 


(3.40) 

(3.41) 


P  ( e  I  or  re  1 )  =  0.5erfc 


/(2Ed/n) 


Au)T 


2n-l 


V— -  |  uw  I  I  I  £.11-1 

2Ei/n  sine  \—  Jcos  j_rrel  +  “ £  AuiT 


+0.5  erfc 


V(2Ed/n) 


/2Ei 

2n-l 

Y(— )sinc(Aa)T/2)cos 

Jrel  *  2  i“T  _ 

(3.42) 
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Since  Yre-|  is  uniformly  distributed  on  [0,2ti]  (due  to  (3.24)  and  the 
theorem  proven  by  Scire  in  [17])  we  can  determine  the  probability  of 
error  by  evaluating 


P(e) 


2n 

/  P(e|rre, 


0 


>  ^rel 


(3.43) 


where  P(e|Yrel)  ‘is  9’ven  by  (3.42)  above. 

Several  simplifications  can  be  made  before  evaluating  (3.43). 

Since  the  integrand,  P(e|Yrel)>  periodic  in  Yrel  with  period  2n  and 
the  integral  is  over  a  full  period  of  Yrel *  both  of  the  two  error 
functions  which  compose  the  integrand  contribute  equally  to  the  value  of 
the  integral.  The  fact  that  the  integration  is  over  a  whole  period  of 
Yrei  allows  us  to  remove  the  factor  of  AiuT/2  from  the  argument  of  the 
cosine  in  the  integrand  without  affecting  the  result.  With  the  above 
simplifications  and  (3.42)  and  (3.43),  the  expression  for  P(e)  becomes, 

1  2* 

P(e)  =  (jp  /  erfc  W(Kd7nJ 

0 

t  AtuT 

+  /  (2E./n)  sine  (— )  cos(Yre1)>  dYpe1  .  (3.44) 


We  define  the  signal  energy-to-noise  density  ratio  at  the  detector  input 
by 

SNR  =  Ed/o,  (3.45) 

and  similarly  define  the  INR  at  the  detector  input  by 

INR  =  Ei /n.  (3.46) 
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These  definitions  can  be  used  together  with  (3.44)  to  calculate 
P(e)  for  various  values  of  SNR  and  INR  by  numerically  evaluating  the 


required  inteqral. 

Figure  3.4  contains  several  curves  which  illustrate  the  detector 

performance  when  the  interference  frequency  offset  is  zero  (i.e.  AwT=0). 

Each  curve  is  a  plot  of  P(e)  vs.  SNR  for  a  fixed  value  of  signal -to- 

★ 

interference  ratio,  SIR  (SIR=Ed/E. )  .  Note  that,  although  this  figure 
is  very  similar  to  that  shown  by  Rosenbaum  [10],  the  approach  taken  in 
the  derivation  is  different  here  as  discussed  previously.  From  this 
figure  we  see  that,  for  a  fixed  SIR,  P(e)  decreases  as  SNR  increases. 
Similarly,  for  a  fixed  SNR  value,  P(e)  increases  as  SIR  decreases.  We 
find  very  poor  performance  for  SIR  values  of  0  dB  and  less. 

Figure  3.5  shows  the  same  data  which  are  plotted  in  Figure  3.4; 
however,  in  Figure  3.5  the  parameterization  has  been  changed.  In  Figure 
3.5  P(e)  is  plotted  vs.  INR  for  several  fixed  values  of  SNR.  Although 
Figures  3.4  and  3.5  show  the  same  data,  one  or  the  other  is  often  more 
useful  when  trying  to  understand  a  particular  facet  of  the  system 
performance.  In  the  remainder  of  this  report  we  shall  often  present 
results  in  similar  pairs  of  plots,  one  showing  P(e)  vs.  SNR  (for  several 
fixed  SIRs)  and  one  showing  P(e)  vs.  INR  (for  several  fixed  SNRs). 


♦The  value  of  SIR  is  constant  for  each  curve  in  Figure  3.4.  Therefore, 
it  is  assumed  that,  as  we  move  along  a  curve  by  varying  the  SNR,  the 
interference  power  simultaneously  varies  in  order  to  keep  the  SIR 
fixed. 
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It  is  interesting  to  examine  P(e)  for  the  ideal  BPSK  detector  as  a 
function  of  SINR  as  defined  by  (2.31).  Since  the  denominator  in  this 
expression  for  SINR  contains  two  terms,  it  is  clear  that  different 
combinations  of  interference  and  noise  power  can  lead  to  the  same  value 
of  SINR.  However,  from  (3.44)  we  see  that  the  ideal  detector  is  not 
equally  susceptible  to  interference  power  and  total  input  noise  power. 
From  this  equation  we  see  that  the  detector  performance  is  dependent 
upon  the  noise  PSD,  the  interference  energy-per-bit-interval  ,  and  the 
interference  frequency.  Therefore  we  might  suspect  that  two  different 
combinations  of  the  input  signals  which  have  equal  SINR  values  might 
produce  different  P(e)  values.  Figure  3.6  shows  that  this  can  indeed 
happen.  This  figure  shows  how  P(e)  varies  with  SINR  for  two  different 
SNR  values.  The  interference  power  was  allowed  to  vary  for  each  curve 
but  the  noise  bandwidth,  noise  PSD  level,  and  the  interference  frequency 
were  held  fixed  for  each  curve.  From  this  figure  we  see  that  vastly 
different  error  rates  can  occur  for  the  same  SINR  value.  For  example, 
for  a  SINR  value  of  -7  dR,  there  is  a  difference  of  7  orders  of 
magnitude  between  the  two  P(e)  curves. 

In  general  we  find  that  P(e)  is  a  function  of  desired  signal  power, 
interference  power,  interference  frequency,  noise  PSD,  and  noise 
bandwidth.  If  all  but  one  of  these  parameters  is  held  constant  and  the 
remaining  parameter  is  allowed  to  vary,  we  find  that  P(e)  does  decrease 
monotonically  as  SINR  increases.  However,  in  the  general  case,  we  must 
calculate  P(e)  for  a  given  scenario  using  (3.44).  The  inability  to 
accurately  predict  system  performance  using  SINR  alone  was  the  prime 
motivation  for  this  study. 


P (E)  VS.  SINR 


SINR  COB) 


Figure  3.6.  Performance  of  ideal  BPSK  detector  vs.  SINR  for  2  SNR 
values  (awT  =  5,  noise  bandwidth  =  10/T,  INR  varying). 
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0.  PERFORMANCE  CALCULATIONS  FOR  COMBINED  LHS 
ARRAY /BPSK  DETECTOR 

In  this  section  we  combine  the  results  of  the  previous  two  sections 
to  determine  the  overall  performance  of  the  BPSK  communication  system 
which  uses  the  LMS  array.  We  first  present  a  set  of  variables  (such  as 
SNR,  signal  arrival  angles,  etc.)  which  completely  specify  the  problem. 
We  then  calculate  the  adaptive  array  performance  using  the  derivations 
in  part  B  of  this  chapter.  Once  we  have  determined  the  power  of  each 
signal  at  the  array  output  (or  PSD  in  the  case  of  the  noise)  we 
calculate  P(e)  of  the  system  using  the  results  of  the  previous  section. 
Several  scenarios  are  then  examined  in  order  to  obtain  both  a 
qualitative  and  quantitative  understanding  of  the  system  performance. 

We  first  present  a  set  of  six  variables  which  specify  the  problem. 
These  variables  are: 


(Ed)in/No  =  t,ie  SNR*  at  each  element  input. 

Efj  is  the  desired  signal  energy -per-bit  at  each  element 
input  and  N0/2  is  the  2-sided  noise  PSD  (in  W/Hz). 

(Ei)in/No  =  at  each  element  input. 

E-j  is  the  interference  signal  energy-per-bit-interval 
at  each  element  input  and  N0/2  is  as  defined  above. 

9(j  =  the  desired  signal  arrival  angle  as  defined  in  Section 

2. 


9-j  =  the  interference  arrival  angle  as  defined  in  Section  2. 


*In  this  report  we  define  the  SNR  and  INR  as  (Ed)in/N0  and  (Ei)in/N0, 
respectively.  These  ratios  are  actually  energy-per-bit-interval  to 
noise  PSD  ratios.  In  most  previous  adaptive  array  work  SNR  and  INR  are 
defined  as  Ad2 /a2  and  Ai2/o2.  We  use  the  former  definition  since  these 
are  the  forms  that  SNR  and  INR  typically  appear  in  the  expressions  for 
the  detector  P(e). 
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k  =  the  array  bandwidth  factor. 

k=BA/(2/T)  where  is  the  input  RPF  bandwidth  (in  Hz.) 
and  T  is  the  bit  duration.  Thus  k  is  the  ratio  of  the 
input  noise  bandwidth  of  the  array  to  the  width  of  the 
main  peak  in  the  desired  signal  PSD.  Figure  3.7 
illustrates  this  definition. 

Au)T  =  the  change  in  phase  of  the  interference  signal  with 

respect  to  the  desired  signal  during  each  bit  interval 
(see  section  C  of  this  chapter). 


To  calculate  the  array  weights  we  only  need  to  specify  the  ratios 
Ad2  Ai2 

r  and  ~T~.  These  ratios  can  be  calculated  given  the  first  five  of  the 
o  o 

variables  listed  above.  We  first  note  that  the  noise  power  at  each  RPF 


output  is  given  by 


a2  =  kN0(2/T).  (3.47) 

We  solve  this  expression  for  N0, 
a2T 

No  =  IF  (3.4ft) 

We  now  write  the  desired  signal  energy  at  each  element  input  as, 

(-d)in  =  No(Ed)in/No*  (3.49) 

Substitution  of  (3.48)  into  (3.49)  yields 
a2T 

(^d)in  =  (  2Y  ( -d)in/No  (3.SD) 

which  we  can  rearrange  as 
( -d)in  j-  (^d)in 

To2  =  2k  NQ  (3. SI) 
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Since  the  desired  signal  power  (Pj)  at  each  BPF  output  is  equal  to 
(Ed)in/T,  (3.51)  becomes 

(pd)in  1  (^d)in 

a2  =  2k  NQ  ’  (3.52) 

or  equivalently, 

V  1  (Ed)in 

a2  =  2k  NQ  •  O-53) 

A  similar  derivation  shows  that. 


Ai2  1  (^i)in 


Aq2  A,2 

Having  determined  and  we  can  apply  the  methods  discussed 

in  section  B  of  this  chapter  to  determine  the  desired  signal  and 
interference  powers  at  the  array  output.  The  noise  PSD  at  the  array 
output  is  given  by  (3.14). 

Once  we  have  determined  the  signals  at  the  array  output  we  can  use 
the  results  of  the  previous  section  to  determine  P(e)  for  the  RPSK 
detector  which  follows  the  adaptive  array. 

The  procedure  outlined  above  was  performed  for  several  different 
scenarios.  Figures  3.8  and  3.9  show  typical  results.  These  figures 
show  how  P(e)  varies  with  SNR  and  iNR  for  the  3-element  array  when  the 
desired  signal  arrives  from  the  broadside  direction  (i.e.,  9^=0°)  and 
the  interference  arrival  angle  9  iS  10°.  The  array  bandwidth  factor 
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Figure  3.7.  Ideal  integrator  and  BPF  transfer  functions  (only  positive 
frequencies  shown). 

k,  was  set  equal  to  5  in  each  of  the  cases  shown  in  these  figures  and 
AojT  was  set  equal  to  zero. 

The  array  behavior  is  most  easily  explained  by  examination  of 
Figure  3.9,  For  each  given  SNR  curve  we  see  that,  as  INR  increases, 

P(e)  first  increases  and  then  decreases.  The  maximum  in  the  P(e)  curve 
occurs  typically  at  INR  values  between  10  and  20  dB.  It  is  in  this  INR 
range  where  the  interference  and  noise  powers  are  roughly  equal  at  the 
element  inputs.  At  INR  values  below  approximately  10  dR  the  thermal 
noise  is  dominant  at  the  array  input  and  the  array  pattern  is  relatively 
unaffected  by  the  interference.  Thus  as  INR  increases,  the  interference 
power  at  the  array  output  and  P(e)  both  increase.  As  the  INR  increases 


above  approximately  10  d8  the  array  begins  to  null  the  jammer.  As  the 
INR  is  increased  further  this  null  becomes  deeper  and  a  point  is  reached 
where  any  additional  increase  in  input  interference  power  produces  a 
corresponding  decrease  in  the  interference  power  at  the  array  output. 
Therefore,  at  very  high  INR  values,  a  very  deep  null  is  formed  in  the 
interfering  signal  direction.  At  these  high  INR  values  the  interference 
power  becomes  negligible  at  the  array  output  (when  compared  with  the 
thermal  noise)  and  the  P(e)  curves  each  asymptotically  approach  a 
constant  value. 

It  is  interesting  to  note  that  the  P(e)  value  which  a  particular 
curve  approaches  at  high  INR  values  is  larger  than  the  value  which  this 
curve  has  at  very  low  INR  values.  This  phenomenon  is  easily  explained. 
At  very  low  INR  values  the  array  is  not  affected  by  the  interference  and 
it  can  adjust  its  weights  in  order  to  maximize  gain  in  the  desired 
signal  direction.  At  very  high  INR  values  the  array  behavior  is 
constrained  so  that  it  must  keep  a  deep  null  in  the  direction  of  the 
i nterference. 

Since  the  desired  and  interfering  signals  are  spatially  close,  the 
desired  signal  lies  very  close  to  the  null.  The  3-element  array  does 
not  have  sufficient  resolution  to  simultaneously  steer  a  null  and  a 
pattern  maximum  only  a  few  degrees  apart.  Therefore,  at  high  INR 
values,  even  though  the  interference  power  at  the  array  output  is 
reduced  to  a  negligible  level,  the  desired  signal  power  is  also  slightly 
reduced.  The  resulting  decrease  in  the  SNR  in  the  array  output  produces 
the  larger  value  of  P(e)  at  very  high  INR  values  compared  to  P(e)  at  low 
INR  levels. 
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The  performance  at  high  INR  values  can  generally  be  improved  by 
adding  elements  (and  thus  degrees  of  freedom)  to  the  array  or  by 
increasing  the  spatial  separation  between  the  desired  and  interfering 
signals.  To  determine  the  effects  of  signal  arrival  angle  on 
performance  we  have  calculated  the  expected  system  P(e)  for  interference 
arrival  angles  of  20°  and  80°.  The  results  for  these  cases  are  shown  in 
Figures  3.10  through  3.13.  As  in  the  previous  cases  the  array  bandwidth 
factor  k,  was  set  equal  to  at  5.  Figures  3.10  and  3.11  show  the  results 
for  9^=20°  and  Figures  3.12  and  3.13  show  the  results  for  0j=8O°. 

Several  interesting  results  can  be  seen  from  these  figures.  Prom  a 
comparison  of  Figures  3.9,  3.11,  and  3.13,  we  observe  the  effects  of  3^ 
on  system  performance.  First  we  see  that,  at  very  low  INR  values,  the 
system  performance  is  not  very  sensitive  to  changes  in  9-j.  In  these 
cases  the  effects  of  the  interference  are  negligible  and  the  array  can 
easily  form  a  pattern  maximum  in  the  desired  signal  direction.  At 
moderate  and  high  INR  values  the  performance  improves  (as  expected)  as 
9j  increases. 

From  Figures  3.8  through  3.13  we  see  that  the  interference  is  most 
effective  when  1  ® i -9<j I  small  and  when  the  INR  is  in  the  10-20  dB 
range. 

In  order  to  illustrate  the  effects  of  the  noise  bandwidth  at  the 
element  inputs,  the  calculations  described  above  were  repeated  with  k 
increased  from  5  to  10.  The  results  of  these  calculations  are  shown  in 
Figures  3,14  through  3.19.  The  results  in  these  plots  correspond 
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directly  to  the  data  shown  in  Figures  3.8  through  3.13  (the  only 
exception  being  the  increase  in  k  from  5  to  10). 

From  these  new  plots  we  see  that  the  array  behaves  similarly  to  the 
previous  cases.  However,  when  the  input  bandwidth  is  increased,  the 
height  of  the  peak  in  the  P(e)  vs.  INR  curves  is  increased  and  this  peak 
is  shifted  slightly  to  the  right.  This  behavior  can  he  explained  as 
follows.  With  the  doubling  of  the  input  bandwidth  twice  as  much  thermal 
noise  is  present  at  each  element  input.  As  described  above,  the  array 
does  not  begin  to  null  the  interference  effectively  until  (Pj)-jn  is 
roughly  equal  to  the  input  thermal  noise  power.  Since  the  increased 
bandwidth  results  in  more  noise  power  at  the  array  input,  a  higher  INR 
is  required  before  the  array  begins  to  reject  the  interference.  Hence 
the  slight  rightward  shift  in  the  curves.  Since  the  maximum  output 
interference  power  level  is  higher  in  the  wider  bandwidth  case,  the 
maximum  0(e)  is  also  higher  in  this  case.  In  light  of  these  results  we 
conclude  that  the  system  performance  is  significantly  degraded  when  we 
doubled  the  input  bandwidth. 

We  have  now  shown  how  the  BPSK/LMS  array  system  is  affected  by  SNR, 
INR,  signal  arrival  angles,  and  array  bandwidth.  Next  we  examine  how 
the  interference  frequency  affects  the  system.  The  contribution  of  the 
interference  to  the  integrator  output  (i.e.  X-j(n))  is  given  by  Equation 
(3. .28).  From  this  equation  we  see  that  Xn(n)  will  be  zero  for  those 
frequencies  at  which  si nc(A<wT/2)  is  zero.  At  these  points  (i.e.,  when 
AajT*2nir  where  n  is  any  non-zero  integer)  the  detector  performance  will 
not  be  degraded  by  the  interference. 
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Figures  3. 2D  and  3.21  show  the  system  performance  when  AusT=2tt.  In 
these  figures  the  array  bandwidth  factor  k  was  set  equal  to  10  and  the 
interference  arrival  angle  0^  was  set  equal  to  10°.  Since  we 
established  in  section  C  of  this  chapter  that  the  detector  is  not 
affected  by  the  interference  with  AtaT=2ir  we  might  be  surprised  to  see 
that  Figures  3.20  and  3.21  show  that  the  system  performance  becomes 
worse  as  the  interference  power  increases.  This  behavior  should  not  be 
surprising,  however,  since  the  adaptive  array  which  precedes  the 
detector  must  change  its  pattern  to  null  the  interference  at  high  INR 
levels.  As  the  array  nulls  this  interference,  the  SNR  at  the  array 
output  decreases  due  to  the  constraints  placed  on  the  array  by  the 
nulling  requirement  (this  phenomenon  is  discussed  above).  The  decreased 
performance  at  high  INR  levels  is  caused  by  this  decrease  in  SNR  at  the 
array  output. 

It  is  interesting  to  note,  however,  that  there  are  no  relative 
maxima  in  the  P(e)  curves  in  Figure  3.21.  These  peaks,  which  we 
observed  in  previous  P(e)  vs.  INR  curves,  were  caused  by  the 
interference  signal  at  the  detector  input  (i.e.,  the  array  output)  at 
moderate  INR  levels.  For  the  case  shown  in  Figure  3.21  the  detector  is 
not  affected  by  the  interference  and  therefore  there  are  no  peaks  in  the 
curves.  In  these  special  cases  (where  /wT=2nir)  SNR  and  not  SINR  is  the 
most  meaningful  measure  of  system  performance  (since  the  detector  is 
immune  to  the  interference).  We  also  find  that,  for  array  bandwidths  of 
interest,  the  SNR  as  defined  by  Equation  (3.51)  is  not  significantly 
dependent  on  the  array  input  bandwidth. 
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Figure  3.10.  BPSK  P(e)  vs.  SNR  for  3-element  array  (9-j=20°,  k=5, 
AojT=0)  . 
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Figure  3.12.  BPSK  P(e)  vs.  SNR  for  3-element  array  (04=80°, 
Au>T=0). 
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P(E)  VS.  SNR 


Figure  3.14.  BPSK  P(e)  vs.  SNR  for  3-element  array  (9^*10°,  k=in, 
AojT=0)  . 
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Figure  3.15.  BPSK  P(e)  vs.  INR  for  3-element  array  (9i=10°,  k=10,  M 

Aa)T=0). 
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Figure  3.18.  BPSK  P(e)  vs.  SNR  for  3-element  array  (9j=80°,  k=10, 
AwT=0) . 
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Figure  3.19.  8PSK  P(e)  vs.  INR  for  3-element  array  (9i=80°,  k  =  lf),  J[ 

AujT=0)  . 
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Figure  3.20.  8PSK  P(e)  vs.  SNR  for  3-eTement  array  (^=10°,  k =10 
AujT=2it  ). 
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If  AujT/2»0  the  sine  function  in  the  expression  for  the  detector 
P(e)  will  be  very  close  to  zero  even  if  AuT*2mr.  Therefore,  for  this 
case,  the  detector  performance  is  not  degraded  by  the  interference. 
However,  the  SNR  at  the  array  output  is  still  reduced  at  high  INR 
levels.  Therefore,  even  though  the  interference  is  outside  of  the 
detector  bandwidth,  the  overall  system  performance  is  still  adversely 
affected  by  the  interference.  The  system  designer  can  limit  this 
susceptibility  to  out-of-band  interference  by  keeping  the  array  input 
bandwidth  as  small  as  possible. 
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CHAPTER  IV 


PERFORMANCE  OF  A  QPSK  SYSTEM  WITH  CW  INTERFERENCE 

A.  INTRODUCTION 

In  this  chapter  we  examine  the  effects  of  CW  interference  on  a 
quadrature  phase  shift  keyed  (QPSK)  communication  system  that  uses  an 
LMS  array.  The  organization  and  analytical  techniques  of  this  chapter 
are  very  similar  to  those  of  the  previous  chapter.  The  equations 
describing  the  QPSK  and  BPSK  systems  are  so  similar  that  we  shall  often 
use  the  results  from  derivations  in  Chapter  III  with  only  minor 
modifications. 

Figure  4.1  shows  a  typical  QPSK  waveform.  Each  of  the  four  possible 
QPSK  symbols  represents  two  bits  of  information.  We  shall  be  careful  to 
use  the  word  "symbol"  and  not  "bit"  to  describe  a  burst  of  energy  at  one 
of  the  signalling  phases.  In  Chapter  III  we  could  use  the  terms  "bit" 
and  "symbol"  interchangeably  since,  for  BPSK  signalling,  each  symbol 
represents  only  one  bit  of  information. 

A  method  for  generating  a  reference  signal  for  an  LMS  array  using 
QPSK  signalling  is  described  by  Winters  in  [31.  This  reference- 
generation  technique,  like  that  discussed  in  Chapter  III  for  BPSK, 
relies  on  spread-spectrum  coding  of  the  desired  signal  in  order  to 
derive  a  reference  signal  from  the  array  output.  We  assume  that  a 
perfect  reference  signal  (one  which  is  identical  to  the  incoming  desired 
signal)  is  available  to  the  array. 
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Figure  4.1.  A  typical  QPSK  waveform  with  bit-pair  assignments. 

B.  ARRAY  PERFORMANCE 

The  array  performance  with  QPSK  signalling  is  very  similar  to  that 
for  the  BPSK  case.  The  desired  signal  at  the  input  to  the  first  array 
element  during  the  nth  symbol  interval  is  given  by, 

d(t)  =  Adexp[u>dt  +  4»n(t )  +  *d],  (4.1) 

where  each  variable  is  identical  to  the  corresponding  variable  in 

Equation  (3.1)  except  that,  for  QPSK,  $  (t)  is  equally  likely  to  be  any 

,  it  3ir, 

member  of  the  set  { 0 ,  ~,  it,  ~}  during  each  symbol  interval. 

It  is  straightforward  to  show  that  the  autocorrelation  function 
(and  thus  the  PSD)  of  the  QPSK  signal  is  the  same  as  that  for  BPSK.  We 
use  the  same  interference  and  noise  models  that  were  used  for  *he  BPSK 
system  in  Chapter  III.  Therefore  the  covariance  matrix  t,  the 
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reference  correlation  vector  S,  and  thus  the  array  weights,  are 
identical  to  those  derived  for  BPSK  in  Chapter  III.  The  desired  signal, 
interfering  signal  and  the  noise  PSD  at  the  array  output  are  given  by 
(3.12),  (3.13),  and  (3.14)  respectively,  the  only  exception  being  that 
<t>n(t)  can  assume  four  values  in  the  QPSK  case. 

Since  the  array  calculations  for  the  BPSK  and  OPSK  systems  are  the 
same,  the  only  difference  in  performance  for  these  two  systems  is  due  to 
the  difference  in  detector  performance.  In  the  next  section  we 
determine  the  performance  of  the  QPSK  detector. 

C.  QPSK  DETECTOR  PERFORMANCE  WITH  CW  AND  AWGN  INTERFERENCE 

In  this  section  we  derive  the  performance  of  an  ideal  QPSK  detector 
when  the  QPSK  desired  signal  is  corrupted  by  CW  interference  and  AWGN. 
The  QPSK  detector  shown  in  Figure  4.2  [201  is  ideal  if  the  only 
undesired  signal  is  AWGN.  This  detector  is  essentially  two  BPSK 
detectors  in  which  the  reference  signals  at  the  two  multipliers  are  in 
phase  quadrature.  At  the  end  of  the  n^  symbol  interval  the  integrator 
outputs,  X(n)  and  Y(n),  are  examined  and  a  symbol  decision  is  made  base 
dupon  the  values  at  these  outputs.  The  decision  regions  for  the 
different  symbols  are  shown  in  Figure  4.3. 

In  order  to  determine  the  probability  of  a  symbol  error  we  must 
first  calculate  the  joint  statistics  of  the  integrator  outputs,  X(n)  and 
Y(n).  We  note  that  the  process  by  which  X(n)  is  derived  is  identical  to 
the  process  by  which  X(n)  was  derived  for  the  ideal  BPSK  detector  of 
Figure  3.3.  Therefore,  from  (4.1)  and  (3.37)  the  expression  for  X(n)  is 


fi6 


sin  (a»dt  +yd) 


Figure  4.2. 


The  ideal  QPSK  detector. 


Figure  4.3.  QPSK  decision  space. 
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immediately  found  to  be 


X(n )  =^(  ~J  T  cos[<j>n(t)] 


_  Iy i  _  ,Aufl\  ,  2n-l 

+  ~ 2  I  sinc(— )cos(yrel  +  — —  AoiT)  +  Xfl(n)  (4.?) 

IT  3'T 

where  i>n(t)  is  equally  likely  to  be  any  member  of  the  set  { 0,  it,  — } 

during  any  symbol  interval.  The  noise  term  X  (n)  is  normally 

n 

nT  n 

distributed  with  zero  mean  and  variance  “where  ~  is  the  noise  PSD  at 
the  array  output. 

We  next  examine  Y(n)  which  is  given  by. 


Y(n)  =  Yd(n)  +  Y^n)  +  Yn(n)  (4.3) 

where, 

nT 

Yd(n)  =  /  v' 2Pd  sin(iuc|t  +  Yd)cos[>dt  +  n ( t )  +  Ydl^t,  (4.4) 

(n-l)T 


nT 

Yi(n)  =  /  sin[udt  +  Ydlcoslj^t  +  Yildt,  (4.5) 

(n-l)T 

and , 

nT 

Yn(n)  =  /  sn(t)sin[o)dt  +  YdM.  (4.5) 

(n-l)T 


As  in  Chapter  III  we  evaluate  each  of  these  three  integrals  and 
then  combine  the  results  to  determine  Y(n). 

It  is  easily  shown  using  techniques  very  similar  to  those  used  to 
derive  (3.22)  that  Yd(n)  reduces  to, 


d 

Yd(n)  “  "Vi-  T  sin^n(t^* 


(4.7) 


The  second  integral  to  be  evaluated,  Yi(n),  is  given  by, 

(V\  nT 

Yi(n)  =  V?“  ^  sin[(ud  +  Wi)t  +  yd  +  Y,-] 

(n-l)T 

+  sin[(u>d  -  uji)t  +  Yd  -  Y^dt. 


(4.8) 


After  we  drop  the  double  frequency  term  we  have. 


Pi  nT 


Mn)  V~2  f  sin(Aut  +  Ypel  )dt, 
(n-l)T 


(4.9) 


where  Aw,  Ypel »  ancl  are  defined  as  in  Chapter  III.  Evaluating  this 
integral  we  have 


i  ,1 


Yi  (n)  =  v  ~2  (^)Ccos(nAwT  +  Trei )  -  cos(n-l)AwT  +Ype1)l.  (4.10) 


Standard  trigonometric  identities  can  be  used  to  simplify  this  equation 
to , 


pi 


Yi  (n )  =~V  ~~2  Tsinc  (^-)  sin  (Yfel  +  AwT) 


(4.11) 


We  denote  the  noise  signal  at  the  integrator  output  at  the  end  of 
the  n**1"1  symbol  interval  by  Yn(n).  Using  an  analysis  nearly  identical  to 
that  used  to  derive  the  statistics  of  X  (n)  as  in  Chapter  III  it  is 

straightforward  to  show  that  Y  (n)  is  a  zero  mean  Gaussian  random 

n '  ' 
nT 

variable  with  variance  ~. 

Furthermore,  if  we  assume  that  the  PSO  of  the  noise  at  the  detector 
input  is  symmetric  about  wd,  then  Xn(n)  and  Yn(n)  will  be  independent  of 
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each  other  at  the  end  of  each  symbol  interval.  This  property  is  easily 
shown  using  the  standard  decomposition  of  narrowband  noise  into 
quadrature  components  [22].  The  values  of  Xn(n),  Yn(n),  X(m),  and  Y(m) 
will  be  uncorrelated  (and  thus  independent)  for  m*n. 

The  desired,  interfering,  and  noise  signals  at  the  integrator 
output  can  now  be  combined  using  (4.3),  (4.7)  and  (4.11)  yielding. 


Y(n) 


n)  =  Tsin 


C*n(t)] 


_  F~\  ,Au>T.  .  2n-l 

h  ,  ~2  T  sine  (—  )sin(Ype1  +  ~  AwT)  +  Yn(n) 


(4.13) 


where  Yn(n)  is  a  zero  mean  normal  random  variable  with  variance  ~. 

The  probability  of  a  symbol  error  can  be  calculated  using  the 
expressions  for  X(n)  and  Y(n)  given  by  (4.2)  and  (4.13),  at  1  the 
decision  regions  shown  in  Figure  4.3.  If  we  assume  equally  likely 
symbols,  then  P(e)  is  given  by 

P(e)  =  |  [P[e|<j>n(t)=0]  +  P[e|<t>n(t)  =  j] 

+  p[e|  >n(t)=Ti ]  +  P[e 1 4>n (t )  =  ”]]  (4.14) 

The  symmetry  of  the  decision  regions  and  the  distribution  of  X(n) 
and  Y(n)  insure  that  the  probability  of  a  symbol  error  is  independent  of 
which  symbol  is  transmitted  (for  equally  likely  symbols).  Therefore  the 
probability  of  error  is, 

P(e)  =  P[e|*n(t)  =  0]  (4. IF) 

If  t>n(t)  =  0  then  X(n)  and  Y(n)  become. 
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X  (n ) 


,  AwT.  . 

Tsinc(— )cos(yre1 


+ 


2n-l 

2 


A<*>T)  +  Xn(n),  (4.16) 


and , 


,  ,  /pi  ,  Au)T 

Y  (n)  ="^  —  Tsinc(— “)sin(Y 


rel 


2n-l 

+  —  AcuTJ 


+  V">- 


(4.17) 


The  problem  of  calculating  P(e)  for  QPSK  given  X(n)  and  Y(n) 
very  similar  to  a  problem  addressed  by  Rosenbaum  [18]  and  we  shall 
similar  notation  and  procedures  for  ease  of  comparison.  We  first 


nq 

normalize  X(n)  and  Y(n)  by  dividing  each  by*^'—  T  yielding 


is 

use 


,  2n-l 

X  =  1  +  bcos(yrel  +  ~2  AwT)  +  nx, 

and , 

2n-l 

Y'  =  bsin(yrel  +  -J—  AiuT)  +  ny, 
where, 

f*T  .  AwT 

b  -pj  Sine  (— ) , 


(4.18) 


(4.19) 


(4.20) 


and  nx  and  ny  are  independent  and  identically  distributed  normal  random 

variables  with  zero  mean  and  variance, 

2  (nT/4) 

°xy  =  (PdT2/2) 

=  (2Ed/n ) *  (4*?n) 


The  normalized  variables  V  and  Y'  are  shown  on  the  decision  space 
diagram  in  Figure  4.4. 

Since  V  and  Y'  are  independent  jointly  Gaussian  random  variables, 
their  joint  probability  density  function  (PDF)  conditioned  on  Yre]  is. 
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Figure  4.4.  Normalized  decision  space  showing  composite  received 
si gnal . 


Vr-r|Yrel>  ”i^r  «P 


xy 


-TJ-  [(X'-bsiny)2  +  (Y ' -1-bcosy  )2  ] 
xy 


where , 


.AuT 


Y  =  Y 


rel  ;  2 


(“T). 


(4.21) 


(4.2?) 


After  eliminating  the  dependence  of  this  PDF  on  y  by  integrating 
the  product  of  f  (X' »Y' I Y rei )  and  ~  (the  PDF  of  the  uniformly 
distributed  variable  y)  we  have. 
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f(x'.Y')  -  (d-)2«P  [ir  C(X'2  *  (X'-D2  +  b23] 

Ay  ^  xy  xy 

2n  b  _ _ 

*  /  exp  [—  [/X'2  +  (Y-l)2  cos(Y+r)l]dY,  (4.23) 

o  xy 


-1  (Y'-l) 

where  r=-cot  ~t  is  not  a  function  of  y. 

This  expression  can  be  simplified  after  we  recognize  the  integral 
to  give  us 

f -  (wt-)  «p  iir  w2  *  (y,-i>2  ♦  i>2i] 


xy’ 


xy 


'xy 


*  In  {A“/X'2  +  (Y'-l)2  },  (4.24) 

xy 

where  I0  is  the  modified  Bessel  function  of  the  first  kind  of  order 
zero. 


We  can  now  change  from  rectangular  to  polar  coordinates  using  the 
transformation , 

X  =  r  sin(ot),  ( 4 . 25 ) 

and , 

Y  =  rcos(a).  (4.2ft) 

The  geometrical  interpretation  of  r  and  a  is  shown  in  Figure  4.4. 

To  implement  the  transformation  we  perform  the  indicated  change  of 
variables,  and  multiply  the  resulting  function  by  the  Jacobian  of  the 
transformation,  r.  We  can  then  eliminate  the  dependency  on  r  by 
integrating  the  density  over  all  possible  values  of  r.  The  resulting 
expression  will  be  the  probability  density  function  of  the  angle  a. 

This  PDF  is  given  by, 
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00 


-1 


[r2  +  b2  +  1  -  2rcosa]j 


f 

a 


(a) 


1 


/  exp  [ 
o 


I  {  ~^z~  /r2  +  1  -  2rcosa}  rdr. 


xy 


(4.27) 


The  probability  of  error,  P(e)  will  be  given  by  the  probability 

-IT  TT 

that  a  does  not  lie  between  “and  “when  =  0.  Thus  P(e)  is  given 

by, 


P(e)  =  f 


-IT /4  TT 

J  fa(a)d«  +  I 
-IT  ir/4 


(4.28) 


From  the  symmetry  of  the  integrand  it  is  seen  that  (4.28)  can  be 
rewritten. 


p(e)  =  2  /  fa(a)da  (4.29) 

tt/4 

Equation  (4.29)  was  evaluated  numerically  for  several  SNR  and  INR 
values.  Figures  4.5  and  4.6  show  results  of  these  calculations  for 
AwT  =  0.  Simpson's  rule  was  used  to  evaluate  the  integral  where  the 
Bessel  function  was  evaluated  by  a  combination  of  polynomial  and 
asymptotic  approximations.  From  Figures  4.5  and  4.6  we  see  that  these 
curves  are  similar  to  the  curves  shown  in  Chapter  III  for  BPSK. 

However,  P(e),  the  probability  of  a  symbol  error  is  higher  for  OPSK.  We 
should  remember  that  each  symbol  in  the  QPSK  symbol  stream  represents 
two  bits  of  information.  Thus,  for  QPSK,  it  is  possible  for  one  symbol 
error  to  produce  two  bit  errors.* 


*  The  probability  of  two  bits  being  in  error  for  a  given  symbol  is  very 
small  if  the  bit  pair  assignments  for  the  decision  regions  are  chosen 
properly.  If  they  are  chosen  so  that  only  one  bit  differs  for  any  two 
adjacent  regions  then  both  X(n)  and  Y(n)  must  be  simultaneously  large 
to  produce  a  double  bit  error. 
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It  should  be  noted  that  the  detector  shown  in  Figure  4.2  is  the 
ideal  detector  for  M-ary  PSK  with  AW (IN.  The  only  change  required  is  the 
division  of  the  decision  space  into  M  regions.  The  probability  of  error 
for  M-ary  PSK  is  given  by  (4.29)  with  the  lower  limit  of  integration 

IT 

changed  to  ~. 

D.  PERFORMANCE  OF  THE  COMBINED  LMS  ARRAY/QPSK  DETECTOR 

In  this  section  we  combine  the  results  of  the  previous  two  sections 
in  order  to  determine  the  overall  performance  of  a  QPSK  detector  that  is 
preceded  by  an  LMS  adaptive  array.  Again  in  this  section  we  find  that 
the  required  derivations  closely  parallel  those  of  the  correspondi ng 
Section  in  Chapter  III.  In  fact,  the  results  from  the  beginning  of 
section  D  of  Chapter  III  apply  exactly  through  Equation  (3.56).  Thus 
the  adaptive  array  calculations  are  identical  for  BPSK  and  QPSK 
signalling.  The  only  difference  in  the  performance  of  the  two  systems 
is  caused  by  the  differences  in  the  BPSK  and  QPSK  detectors. 

The  performance  of  the  LMS  array/OPSK  detector  was  calculated  for 
several  scenarios.  Results  are  shown  in  Figure  4.7  through  4.2(1.  The 
scenarios  examined  correspond  one-to-one  with  the  curves  shown  in 
Figures  3.7  through  3.20.  Figures  4.7  and  4.8  show  the  system 
performance  with  9<j  =  0°,  9j  =  10°,  the  array  bandwidth  factor  k  =  5, 
and  the  interference  phase-off set-per-symbol  (AcoT)  equal  to  zero. 

Figures  4.9  through  4.12  show  similar  results  for  9-j  values  of  20 
and  80  degrees.  From  these  figures  we  see  that  the  QPSK  performance 
curves  are  shaped  similarly  to  those  for  BPSK.  However,  the  QPSK  symbol 
error  probability  is  higher  in  each  case. 
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(Old)  001 


Figure  4.7.  QPSK  P(e)  vs.  SNR  for  3-element  array  (9 ^  =10° ,  k=5,  AjT=0). 


I 

I 


P (E)  VS.  INR 


Figure  4.10.  QPSK  p(e)  vs.  INR  for  3-element  array  (^=20°,  k=5, 
Aa)T =0)  . 
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LOG  IP  IE)  J 


P  (E)  VS.  SNR 


SIR  >  OdB 


-15,-30 


6.  8.  10.  12.  11.  16.  18.  i 

SNR  (OB) 


Figure  4.11.  QPSK  P(e)  vs.  SNR  for  3-element  array  (•>1=80 


Figures  4.13  through  4.18  show  the  system  performance  with  the 
array  bandwidth  factor  k,  increased  to  10  with  all  other  variables  the 
same  as  in  Figures  4.7  through  4.12.  Figures  4.19  and  4.20  show  the 
performance  with  the  interference  phase-shi ft-per-symbol  AujT  set  equal 
to  it.  Again  in  these  cases  we  see  that  the  system  performance  is 
similar  to,  but  slightly  worse  than,  the  RPSK  system  of  Chapter  III. 

The  qualitative  descriptions  of  the  system  behavior  are  identical 
for  the  BPSK  and  QPSK  cases  and  the  discussion  at  the  end  of  Chapter  III 
applies  directly  to  the  QPSK  case.  The  only  difference  in  the  two  cases 
is  the  uniformly  poorer  performance  of  the  QPSK  case.  In  return  for 
this  performance  degradation,  QPSK  offers  a  bit  rate  double  that  of  RPSK 
for  a  fixed  signal  bandwidth. 
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7  i - i - 1 - 1 - 1  i  i  r  i 

6.  8.  10.  12.  IN.  16.  18.  20.  22. 

SNR  (DB) 


Figure  4.13.  QPSK  P(e)  vs.  SNR  for  3-element  array  (9^10°,  k=10, 
AioT=0) . 
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LOG  IP  IE) ) 


P  (E)  VS.  SNR 


Figure  4.15.  QPSK  P(e)  vs.  SNR  for  3-element  array  (9^20°,  k=in, 
AojT=0)  . 
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P  (E) 


II.  16.  18.  20.  22 

iNR  (OB) 


3-element  array  (9-j=RD°,  k=10 


LOG  IP  (EM 


P(E)  VS.  SNB 


Figure  4.19.  QPSK  P(e)  vs.  SNR  for  3-element  array  (0^=10°,  k=10, 
AujT=2ti  ). 
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Figure  4.20.  QPSK  P(e)  vs.  INR  for  3-element  array  (9^=10°,  k=10, 
£m)T=2. ) . 
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CHAPTER  V 


PERFORMANCE  OF  A  DPSK  SYSTEM  WITH  CW  INTERFERENCE 

A.  INTRODUCTION 

In  this  chapter  we  examine  the  effects  of  CW  interference  on  a 
(binary)  differential  phase-shi ft -keyed  (DPSK)  communication  system  that 
uses  an  IMS  array.  DPSK  and  BPSK  systems  are  very  similar;  however, 

DPSK  systems  do  not  require  a  phase-coherent  reference  signal  at  the 
detector.  Although  the  transmitted  DPSK  signal  and  PSD  look  identical 
to  those  of  BPSK  (see  Figure  3.1)  the  method  of  modulating  the  carrier 
with  the  data  bits  is  different  for  the  two  systems.  In  BPSK  systems, 
the  phase  of  the  transmitted  waveform  during  each  symbol  interval  is 
determined  directly  by  the  bit  to  be  transmitted,  one  phase  value  for  a 
"0"  bit  and  another  for  a  "1”  bit.  In  DPSK  systems  the  data  is  encoded 
differentially.  That  is,  a  "1"  is  transmitted  by  the  introduction  of  a 
phase  reversal  between  two  successive  symbols  and  a  "0"  is  transmitted 
by  the  absence  of  a  phase  change  between  the  two  symbols. 

The  DPSK  detector  compares  the  phase  of  each  received  symbol  with 
the  phase  of  the  previous  symbol.  If  the  two  phases  are  nearly  the  same 
the  detector  decides  that  a  "0"  was  transmitted.  If  the  two  phases 
differ  by  more  than  90°,  the  detector  decides  that  a  "1"  was 
transmitted.  Of  course  the  DPSK  detector  can  be  used  with  more  than  two 
signalling  phases  by  appropriately  partitioning  the  detector  decision 
space. 


93 


We  first  examine  the  array  performance  with  a  DPSK  desired  signal. 
We  then  calculate  the  detector  performance  and  the  overall  system 
performance. 


8.  ADAPTIVE  ARRAY  PERFORMANCE 


The  binary  DPSK  desired  signal  at  the  output  of  BPF1  is  given  by, 
d(t)  =  Adexp[ojdt  +  +n(t)  +  +d]  (n-l)T<t<nT  (5.1) 


where  all  variables  except  4>n (t )  are  defined  as  in  Equation  (3.1). 
During  the  ntfl  bit  interval  4>n (t )  is  given  by. 


♦n<t>  ’ 


•  *(n-l)ft>  +' 


during  "0"  symbols 


during  "l”  symbols. 


(5.2) 


Note  that  for  DPSK,  as  for  BPSK,  the  (modulo  2n)  value  of  4>n(t)  is 
equally  likely  to  be  0  or  it  during  each  symbol  interval.  To  the  LMS 
array,  the  BPSK  and  DPSK  signals  are  indistinguishable.  Therefore  the 
results  of  Section  B  of  Chapter  III  apply  directly  to  the  DPSK  case. 
These  results  can  be  used  to  calculated  the  signal  powers  at  the  array 
output. 

The  difference  in  performance  for  the  DPSK  and  BPSK  systems  is  due 
only  to  the  difference  in  detector  performance.  In  the  next  section  we 
calculate  the  performance  of  the  ideal  DPSK  detector. 
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C.  DPSK  DETECTOR  PERFORMANCE  WITH  CW  AND  AWGN  INTERFERENCE 


In  this  section  we  determined  the  performance  of  a  DPSK  detector 
when  the  desired  signal  is  corrupted  by  both  AWGN  and  CW  interference. 

We  examine  the  effects  of  signal  powers,  arrival  angles,  and  the 
interference  frequency. 

Figure  5.1  shows  an  ideal  DPSK  detector  [23].  The  incoming  signal 
is  split  and  the  two  resulting  signals  are  multiplied  by  sinfw^t)  and 
cos(wdt).  These  locally  generated  reference  signals  are  not 
phase-coherent  with  the  incoming  desired  signal  (due  to  the  random 
phase  angle  yj  of  the  desired  signal).  The  multiplier  outputs  are 
integrated  over  the  symbol  interval  to  produce  two  random  variables, 

X(n)  and  Y(n).  At  the  end  of  the  n^*1  symbol  interval,  X(n)  and  Y(n)  are 
compared  with  X(n-l)  and  Y(n-l)  and  a  symbol  decision  is  made  based  upon 
thr  value  of  X(n)X(n-l)  +  Y(n)Y(n-l).  The  decision  rule  is, 

X(n)X(n-l)  +  Y(n)Y(n-l)  >  0  =>  decide  nth  symbol  was  "0", 

X(n)X(n-l)  +  Y(n)Y(n-l)  <  0  =>  decide  nth  symbol  was  "1". 

Figure  5.2  shows  the  DPSK  decision  space  and  typical  values  of  X(n-l), 
Y(n-l),  X(n),  and  Y(n).  From  this  figure  we  see  that  the  decision 
boundary  produces  a  useful  gometric  interpretation.  If  the  absolute 
value  of  the  angle  (6)  between  the  vectors  drawn  from  the  origin  to  the 
points  <X(n-l),  Y(n-1)>  and  <X(n),  Y(n)>  is  less  than  90°  a  "0"  decision 
is  made.  If  | 6 |  is  gr.ater  than  90°  a  "1"  decision  is  made. 

We  now  calculate  the  performance  of  the  DPSK  detector  when  the 
detector  input  is  of  the  form. 
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<x(n-l),Y(n-l)> 


Figure  5.2.  DPSK  decision  space. 

s(t)  =  /2PJ  COS(wdt+4>n(t)+Yd)  +/2P7cos(a)i  +Yi)  +  n(t).  (5.3) 

where  all  terms  are  as  defined  in  Chapter  III  except  for  <j>n(t)  which  is 
defined  in  (5.2). 

The  signal  at  the  output  of  the  upper  integrator  in  Figure  5.1  at 
the  end  of  the  nth  symbol  interval  is  given  by, 

nT 

X(n)  =  /  s(t)cos(«dt)  dt 
(n-l)T 

=  Xd(n)  +  X-  (n)  +  Xn(n),  (^-41 
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Figure  5.3.  DPSK  composite  signal, 
where 


nT 

^  cos(u>dt  +  ^(t)  +  yd)cos(o)dt)  dt 

Xd(n)  -  I 

(5.5) 

(n-l)T 

nT 

/2PT  COS (<u . t  +  Yi  )cOS(oidt)  dt 

M")  =  J 

(n-l)T 

(5.6) 

nT 

Xn (n)  "  / 

n(t)  cos(u>dt)  dt 

(5.7) 

(n-l)T 

Similarly  the  output  of  the  lower  integrator  at  the  end  of  the  nth 

symbol  interval  is 
nT 

Y(n)  =  /  s(t)sin(ujdt)  dt 

(n-l)T 

3  V")  +  M")  +  Yn(n),  (5.8) 
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where, 


nT  _ 

Yd(n)  =  /  /2Pd  cos(a)dt  +  *n(t)  +  Yd)sin(a)dt)  dt  (5.9) 

(n-l)T 

nT  _ 

Y. (n)  =  /  /2Pi  cos(u.t  +  )sin(wdt)  dt  (5.10) 

(n-l)T 


nT 

Yn(n)  =  /  n(t)s1n(o»dt)  dt  (5.11) 

(n-l)T 

We  now  determine  Xd(n),  Yd(n),  X^(n),  and  the  distribution  of  the  random 
variables  Xn(n)  and  Yn(n).  We  simplify  Xd(n)  by  first  expanding  the 
integrand  in  Equation  (5.5).  We  then  drop  the  double  frequency  term  and 
obtain  the  expression. 


_  nT 

xd(n)  «  ^Pd/2  /  cos(2u>dt+4>  (t)+yd)  +  cos(<j)  (t)+yd)dt 

(n-l)T 

=  SVJZ  Tcos(*n(t)  +  Yd).  (5.12) 

Similarly  Yd(n)  is  simplified  from  (5.9), 

_  nT 

Yd(n)  =  ^Pd/2  /  s1n(2udt+*  (t)+yd)  -  sin(*  (t)+yd)dt 

(n-l)T 

-  -  Tsin(<j>n(t)  +  Yd).  (5.13) 

The  first  of  the  interference  terms,  X^(n)  is  simplified  from  (5.6)  in 
the  same  manner. 
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nT 

/  COS[(u)d-Ho.  )t+Yi] 

(n-l)T 
+  cos[(«d-«1)t-Yi]  dt 

=  /fTJz  (l/Aio){sin[AuinT-yi]  -  sin(Ato(n-l)T-Yi )}  (5.14) 

where  Aw  is  as  defined  in  the  previous  chapters.  We  simplify  this 
expression  further  using  standard  trigonometric  identities, 

X.  (n)  =  Aj/2  Tsinc(A<uT/2)  cos[(2n-l)(AuT/2)  -  y^.  (5.15) 

A  similar  procedure  shows  that  Y.  is  given  by, 

Y. (n)  =  /pi 7 2  Tsinc  (AuiT/2)  sin[(2n-l) (Aq>T/2)  -  y.].  (5.16) 

We  next  determine  the  distribution  of  Xn(n)  and  Yn(n),  the  noise 
terms.  These  variables  are  the  result  of  processing  that  is  identical 
to  that  used  to  obtain  Xn(n)  and  Yn(n)  for  QPSK  in  Chapter  IV.  Using 
the  results  from  Chapter  IV  we  can  immediately  determine  that  Xn(n)  and 
Yn(n)  are  independent  zero-mean  Gaussian  random  variables  each  with 
variance  nT/4  (where  n/2  is  the  2-sided  noise  PSD  at  the  detector 
input.)  Xn(n)  and  Yn(n)  will  each  be  independent  of  Xn(n+c)  and  Yn(n+c) 
for  all  non-zero  integer  values  of  c. 

Figure  5.3  shows  each  of  the  components  of  X(n)  and  Y(n)  as  well  as 
the  composite  signals.  This  figure  shows  the  decision  regions,  which 
are  separated  by  the  line  at  |S|=90°. 

In  order  to  determine  the  probability  of  error  we  must  determine 
the  probability  that  1 6  j >90°  when  "0"  is  transmitted  and  the  probability 
that  | 5 | <90°  when  "1"  is  transmitted.  Rosenbaum  [18]  addressed  a 


x,  (n)  =>~/2 
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similar  problem;  however,  as  in  the  BPSK  and  QPSK  cases,  Rosenbaum's 
receiver  model  differs  from  ours  and  his  expressions  for  P(e)  must  be 
slightly  modified  for  the  integrate-and-dump  type  detector.  The  model 
discussed  in  [18]  leads  to  a  signal  vector  diagram  very  similar  to  that 
shown  in  Figure  5.3.  The  only  difference  is  that  Rosenbaum  omits  the 
factor  of  sinc(Au)T/2)  for  the  length  of  the  interference  vector.  We 
outline  Rosenbaum's  derivation  of  P(e)  with  the  appropriate 
modifications  of  the  vector  lengths. 

In  order  to  determine  the  P(e)  we  first  define  the  complex  variable 

Z(n), 

Z(n)  =  X(n)  +  jY(n).  (5.17) 

We  assume  that  the  nth  symbol  is  a  "0".  We  then  note  that 
Re[Z(n)Z*(n-l)] 

cos6  =  | Z ( n )  1 1 Z(n-l)  |  *  ^,18) 

Since  the  denominator  is  nonnegative,  the  conditional  probability  of 
error,  given  AwT,  Yd,  and  ,  is, 

P(e|Aa>T,"0",Yd,Yi)  =  P(Re[Z(n)Z*(n-l)]  <  0)  (5.19) 

Since 

Z(n)  +  Z(n-l)  2  |  Z(n)  -  Z(n-l)  2 

Re[Z(n)Z*(n-l)]  =  - - -  -  |  - - -  ,  (5.20) 

we  have, 

P(e|Au)T,"0",Yd,Yi)  =  P( | Z(n)  +  Z(n-l)|  <  |Z(n)  -  Z(n-l) | ) ,  (5.21) 
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We  now  define  the  variables, 
l  =  Z(n)  +  Z(n-l) 

=  CQ(I)  +  Xn(n)  +  Xn(n-1)  +  j[Yn(n)  +  Yn(n-1)]  (5.22) 

and 

8  =  Z(n)  -  Z(n-l) 

-  C0(B)  +  Xn(n)  -  Xn(n-1)  +  J[Yn(n)  -  Yn(n-1)]  (5.23) 

where  C  (Z)  and  CQ(8)  are  the  components  of  z  and  8  not  due  to  the 
Gaussian  noise, 

c0(£)  «  Xd(n)  +  Xd(n-1)  +  Xi  (n)  +  X^n-1) 

+  j[Yd(n)  +  Yd(n-1)  +  Y.(n)  +  Y.(n-l)]  (5.24) 

C0(3)  =  Xd(n)  -  Xd(n-1)  +  X.(n)  -  X^n-1) 

+  j[Yd(n)  -  Yd(n-1)  +  Y^n)  -  Y.(n-l)]  (5.25) 

where  the  subscript  "0"  indicates  that  we  are  assuming  the  nth  symbol  is 
a  zero.  Evaluating  CQ(z)  and  CQ(8)  we  find, 

C0(s)  =  cos(*n) 

Afl)T  A(i)T 

+  /2P,jT  sinc(~^~)  cos[yd  -  ri  +  AuiT(n-l)]  cos(“£“j 

+  k[-/2Pj  sin(*n) 

Maa  AwT  f  AwT.  «  .  . 

+  /2PjT  sinc(-Y“)  sinCYd  -  Yi  +  AwT(n-l)]  cos(-yj]  (5.26) 
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and 


(5.27) 


We  can  show  that  |i|  and  | @ |  are  independent  Rician  random 
variables.  Therefore,  the  probability  given  in  (5.21)  is  the 
probability  that  one  Rician  random  variable  exceeds  another.  This 
probability  is  given  by  [18], 


1  |C0(E)|  | C0(3 ) | 

P(e|A«T,"0",Yd,Yi)  =  7  [i-Q 


|c0(e)|  I c0(e ) | 

^  *  To  )] 


(5.28) 


nT 


where  a  »  and  Q(A,B)  is  Marcum's  Q  function. 


0(A,B)  =  /  t  exp[-(A2+t2 )/2]  Iq(At)  dt, 
B 


(5.29) 


From  (5.26)  and  (5.27)  we  find  that  I CQ (E ) |  and  | CQ (B ) |  can  be 
written. 


2,  AujT.  ?,Ao)T' , 

Pi  Slnc  (— )cos^(— )] 


1/2 


and , 


|Cq(E)|  = 


AojT  AujT 

22^  sinc(~yj  sin  (""£”) 


(5.30) 


(5.31) 


where  =  [4>n  -  Yd  +  Yi  -(n-l)AwT].  Since  Yd  and  y1  are  1J[0,2*],  50  is 
also  U[0,2tt]  (as  shown  in  [17]). 
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We  now  average  the  expression  for  P(e),  given  by  Equation  (5.28), 
over  the  uniformly  distributed  random  variable  £0  to  determine 
P(e|  AtoT,"0") , 


P(ej  Au)T,"0") 


=  W  [1-Q( 


1C0(E)1  |C0(8)| 
2o  ’  2a 


■)  +  0( 


|C0(S)|  |c0(e ) |  , 


o’ 


(5.32) 


A  similar  analysis  performed  under  the  assumption  that  a  "l"  symbol 
was  transmitted  yields. 


1  2tt  |Ci(E)|  |C!(S)| 

pcitoT.-n-^/  [l-ot-ir- . -i  ) 

o 


♦  Q( 


|Ci(b)|  |Ci(s)l_ 

2a  ’  2a  >>  dCl* 


(5.33) 


where. 


l-  y—  ■  AuT  IT  -Au)T 

|C1(Z)|  = /2T[Pd  +  /PjPT  sine  (— )cos(— y-JcosUj) 


and , 


2,Aa)T\  ?,w-A«Tx ,  i  /? 
+  Pi  sine  (— )cosl— g—)]1'* 


r—  ■  A(l)T  IT 

I C x ( 8 ) |  =  /2p7t  sine  (— )  sin  (— — ) 


(5.34) 


(5.35) 


and  is  U[0,2ir]. 

Finally,  the  expression  for  P(e)  is  given  by, 
P(e|AwT)  =  |  P(e|A«oT,"0")  +  \  P(e|Au>T,"l"). 


(5.36) 
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We  evaluated  P(e)  for  several  values  of  SNR  and  INR  (defined  by 
(3.45)  and  (3.46))  numerically.  Figures  5.4  and  5.5  show  typical 
results  for  AwT=0.  These  results  are  very  similar  to  those  shown  for 
AwT=0  in  [18].  Figures  5.6  and  5.7  show  similar  results  for  AwT=ir/2. 
From  Figures  5.6  and  5.7  we  see  that  the  system  performance  decreases 
slightly  as  AwT  increases. 

The  curves  shown  in  Figure  5.6  depict  significantly  better 
performance  than  that  predicted  by  the  results  of  [18]  for  AwT=m/2. 
There  is  a  simple  explanation  for  this  discrepancy.  For  a  given 
interference  signal  vector  length  (in  Figure  5.3)  the  DPSK  detector  is 
most  prone  to  errors  when  AuiT=ir/2  (±nr).  However,  the  length  of  this 
vector  is  proportional  to  sine (AwT/2) ,  the  factor  which  is  not  included 
in  the  analysis  of  [18],  Therefore,  as  AwT  increases  from  zero  to  tt/2, 
two  processes  occur.  First,  the  detector  sensitivity  to  the 
interference  decreases  due  to  the  sine (AwT/2)  term.  Second,  the 
detector  becomes  more  sensitive  to  the  interference  as  AwT  approaches 
ir/2  due  to  the  geometry  of  the  signal  vectors. 

In  order  to  determine  how  these  two  (often  conflicting)  processes 
interact,  we  examine  how  P(e)  varies  as  a  function  of  AwT  for  few 
representative  SNR  and  INR  levels.  Figures  5.8  and  5.9  show  typical 
results.  Figure  5.8  shows  P(e)  vs.  AwT  for  several  INR  values  when  the 
SNR  is  10  dB.  Figure  5.9  shows  similar  data  for  the  same  INR  values 
with  an  SNR  of  13  dB.  From  these  curves  we  can  see  that  the  highest 
P(e)  values  typically  occur  when  AwT  is  approximately  w/2. 
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0.  PERFORMANCE  OF  COMBINED  LMB  ARRAY/DPSK  DETECTOR 


We  now  combine  the  results  from  the  previous  two  sections  in  order 
to  calculate  the  performance  of  the  DPSK  system.  Figures  5.10  through 
5.25  show  typical  results.  Figures  5.10  through  5.15  show  how  P(e) 
varies  with  SNR  and  INR  with  the  array  bandwidth  factor,  k,  set  equal  to 
5.  As  in  the  previous  chapters,  the  figures  show  the  performance  for 
three  different  interference  arrival  angles,  10,  20  and  80°.  Figures 
5.16  through  5.21  show  similar  results  when  the  array  bandwidth  factor 
is  increased  to  10.  From  these  figures  we  see  that  the  DPSK  system 
performance  is  significantly  poorer  than  that  of  the  BPSK  system 
discussed  in  Chapter  III. 

Figures  5.22  and  5.23  show  how  the  performance  varies  when  the 
interference  and  the  desired  signal  frequencies  are  different.  For  each 
of  these  figures  k  is  10,  9-j  =  10°,  and  AuT  is  w/2.  From  the  previous 
section  we  expect  this  value  of  AujT  to  produce  the  worst  performance. 

Figures  5.24  and  5.25  show  the  performance  for  AwT=2ir  (with  k  =  10 
and  0-j  =  10°).  For  this  value  of  AuT,  the  detector  is  not  affected  by 
the  interference.  Therefore,  as  we  would  expect,  we  observe  the  best 
performance  for  this  value  of  AwT. 

We  find  that  the  system  behavior  with  DPSK  signalling  is  similar  to 
that  with  BPSK  or  QPSK  signalling.  We  find  that,  in  general,  the  symbol 
error  probability  for  DPSK  is  between  that  of  BPSK  and  QPSK.  The  shapes 
of  the  curves  for  each  of  the  different  types  of  PSK  signalling  are  very 
similar.  We  find  that,  for  each  type  of  PSK  modulation,  it  is  desirable 
to  have  a  narrow  array  input  bandwidth. 
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In  Chapters  III,  IV,  and  V  we  have  examined  the  performance  of  PSK 
communication  systems.  We  found  qualitatively  similar  performance  for 
each  of  these  systems.  In  the  next  chapter  we  examine  the  performance 
of  an  FSK  system  and  compare  its  performance  with  that  observed  for  the 
PSK  systems. 
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((3)  <11  001 


P (E)  VS.  SNR 


SNR  (OB) 


Figure  5.4.  Performance  of  ideal  OPSK  detector  with  CW  interference 
(a<dT=0)  . 
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PIE)  VS.  SNR 


SNR  (OB) 


Figure  5.6.  Performance  of  ideal  DPSK  detector  with  CW  interference 
(Au>T=ir/2) . 
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LOG  IP  (EH 


P  (E)  VS.  INR 


Figure  5.7.  Performance  of  ideal  DPSK  detector  with  CW  interference 
(AojT=w/2). 


Ill 


P  (E) 


SNR 


Figure  5.10.  OPSK  P(e)  vs.  SNR  for  3-element  array  (9-j=10°,  k=5, 
AuT=0) . 
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T - 1 - i  I  I  1 

?0.  -10.  0.  10.  20.  SO. 

INR  (OBI 


Figure  5.11.  DPSK  P(e)  vs.  INR  for  3-element  array  (9-j=10o,  k=5, 
AujT=0)  . 
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P  (E)  VS.  INR 


Figure  5.13.  OPSK  P(e)  vs.  INR  for  3-element  array  (9-j  =20 
AwT=0). 
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P(E)  VS.  SNR 


Figure  5.14.  OPSK  P(e)  vs.  SNR  for  3-element  array  (9i=80°,  k=5, 
&u>T=0)  . 
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CHAPTER  VI 


PERFORMANCE  OF  AN  FSK  COMMUNICATION 
SYSTEM  WITH  CW  INTERFERENCE 


A.  INTRODUCTION 

In  this  chapter  we  examine  the  performance  of  a  noncoherent  binary 
frequency  shift  keyed  (FSK)  communication  system  that  uses  an  LMS 
adaptive  array  to  suppress  CW  interference.  Figure  6.1  shows  a  typical 
binary  FSK  waveform.  "1"  and  "0"  bits  are  each  represented  by  a  burst 
of  energy  at  one  of  two  signalling  frequencies  (uj  and  u>0 
respectively). 

Figure  6.2  shows  a  simple  non-coherent  FSK  detector.  In  this 
detector,  the  received  signal  is  passed  through  two  bandpass  filters 
(BPFs),  each  of  which  is  tuned  to  one  of  the  signalling  frequencies. 

The  output  of  each  filter  is  envelope  detected.  The  envelope  detector 
outputs  are  sampled  at  the  end  of  each  bit  interval  and  bit  decisions 
are  made  based  upon  a  comparison  of  these  samples. 

Hudson  [6]  proposed  a  method  of  generating  an  FSK  reference  signal 
for  the  LMS  array.  This  method  requires  that  the  transmitted  bit  stream 
be  encoded  into  a  Markov  symbol  stream  prior  to  transmission.  The 
reference  generation  circuit  uses  the  known  Markov  transition 
probabilities  of  the  symbol  stream  to  make  symbol  predictions.  This 
circuit  uses  the  symbol  predictions  to  generate  a  reference  signal  for 
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the  array  that  is  partially  correlated  with  the  desired  signal.  The 
feasibility  of  this  system  was  experimentally  verified  in  [7]. 

In  this  chapter,  as  in  the  previous  chapters,  we  assume  that  a 
reference  signal  is  available  to  the  array  which  is  a  perfect  replica  of 
the  desired  signal.  Although  the  reference  generation  method  discussed 
above  does  not  generate  a  perfect  reference  signal  (due  to  symbol 
prediction  errors)  the  assumption  of  a  perfect  reference  signal  makes 
the  array  weight  calculations  more  tractable. 


AMPLITUDE 


Figure  6.1.  A  typical  FSK  waveform. 
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B.  ARRAY  PERFORMANCE 


In  this  section  we  use  the  results  from  Chapter  II  to  derive  the 
signal  power  levels  at  the  output  of  the  LMS  array  for  the  case  of  an 
FSK  desired  signal.  We  first  describe  the  real  form  of  the  desired 
signal  and  calculate  its  ACF.  We  then  use  this  real  ACF  to  calculate 
the  complex  ACF  (which  we  will  need  in  the  array  calculations). 

The  desired  signal  during  the  bit  interval  is  given  by, 

d(t)  =  Adcos(u)Qt  +  (6.1) 

during  "0"  bits  and 

d(t)  =  Adcos(a)1t  +  fj),  (6.2) 

during  "1"  bits  where 

Ad  =  the  desired  signal  amplitude, 

Wq  =  the  "0"  signalling  frequency, 

Wj  =  the  "1"  signalling  frequency, 

4<q  *  the  phase  angle  of  the  carrier  during  "0"  symbols, 

t|>j  =  the  phase  angle  of  the  carrier  during  "1"  symbols. 

We  assume  that  and  <|^  are  independent  and  identically 
distributed,  each  ~  IJ[0,2*],  and  that  T  is  the  symbol  duration.  The  FSK 
signal  described  above  might  be  generated  by  switching  between  two 
free-running  oscillators,  one  tuned  to  each  of  the  signalling 
frequencies. 
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We  assume  that  the  incoming  symbol  stream  is  symmetric  (i.e.  "0"s 
and  "l"s  are  equally  likely)  and  that  each  symbol  is  statistically 
independent  of  all  others.  Thus  each  symbol  represents  one  bit  of 
information. 

We  assume  that  the  interference  is  a  CW  signal  that  is  exactly 
aligned  in  frequency  (but  not  necessarily  in  phase)  with  one  of  the 
desired  signal  symbols.  If  we  arbitrarily  choose  as  the  interference 
frequency,  we  have 

?(t)  =  Ai  exp(coQt  +  4»i )  (6.3) 

where  A.,  is  the  interference  amplitude  and  ij>.  is  the  interference  phase 
angle  which  is  ~U[0,2tt]  and  independent  of  both  4^  and  <|^. 

We  now  calculate  the  autocorrelation  function  of  the  desired  and 
interfering  signals.  The  analytic  ACF  for  the  interference  is  given  by 
(3.4).  For  the  desired  signal  we  first  determine  the  ACF  for  the  real 
signal  d(t)  and  then  calculate  Rd(x)  using  the  relationship  [15], 

R~(t)  =  2[Rd(t)  +  jRd(x)],  (6.4) 

where  Rd(x)  is  the  ACF  of  3(t),  Rd(x)  is  the  ACF  of  d(t),  and  the  carat 
denotes  the  Hilbert  transform. 

The  ACF  of  the  desired  signal  as  described  above  was  calculated  by 
Bennet  and  Rice  [25]  for  the  case  of  an  infinitely  long  symbol  stream 
under  the  conditions, 

V  * 
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Figure  6.2.  FSK  detector. 
oj^T  *  irM, 

and, 

ui^T  +  uijT  *  nN, 

where  T  is  the  symbol  duration  and  K,  M,  and  N  are  arbitrary  integers. 
These  conditions  are  not  particularly  restrictive  for  real  systems  since 
the  instability  of  real  oscillators  insures  that  the  inequalities  are 
satisfied.* 

Under  the  above-listed  conditions  the  ACF  for  the  real  signal  d(t) 
is  given  by  [25], 


*  We  assume  that  no  special  effort  is  made  at  the  transmitter  to 
synchronize  the  bit  intervals  with  the  carriers. 
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Ad  (T  -  |t|)  .  Ad  . 

Rd(x)  =  ~ 8  T  (cosujgT  +  coster)  +  —  (coswqT  +  costo^x) 

(6.5) 

for  1 t j <T.  This  result  can  also  be  obtained  by  calculating  the  ACF  of 
two  on-off -keyed  (00K)  waveforms,  one  with  carrier  frequency  u)q  and  the 
other  at  frequency  co^.  Addition  of  these  two  ACFs  yields  Equation 
(6.5). 

Substitution  of  R(j(t)  from  (6.5)  into  (6.4)  yields,  Rd, 

Aj  l  . 

d  M 

Rd  =  2  “  (1  -  pp-  )[(cosajf5x  +  cosajp)  +  j(sino)0x  +  sinwp)] 


-  2t 


JV 


+  e 


JV- 


(6.6) 


The  desired  signal  covariance  matrix  is  calculated  by  substitution 
of  (6.6)  into  (2.20).  Before  writing  the  expression  for  4>d  we  make  some 
simplifying  assumptions  similar  to  those  that  we  made  in  previous 
chapters  for  the  PSK  systems.  First,  we  assume  that  the  interelement 
propagation  delay  is  much  shorter  than  the  symbol  duration,  (i.e. 

Td  «  T).  We  also  assume  that  the  signalling  frequencies  are  close 
enough  so  that  the  interelement  phase  shift  for  the  desired  signal  is 
essentially  independent  of  the  transmitted  symbol  (i.e.,  u)prd  «  (DjT^  = 

*d)- 

These  simplifying  assumptions  allow  us  to  approximate  the  desired 
signal  ACF  by 

RJ(t)  -  A^ejVr  (6.7) 
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Under  the  above  assumptions  the  desired  signal  covariance  matrix 
becomes , 


(6.8) 


The  simplifications  used  to  obtain  the  approximation  (6.8)  do  not 
seriously  affect  the  array  performance  calculations.  Although  these 
simplifications  essentially  cause  us  to  neglect  the  desired  signal 
bandwidth,  it  is  shown  in  [16]  that  the  array  behavior  is  not  strongly 
dependent  on  the  desired  signal  bandwidth  (for  the  bandwidths  of 
interest  here).  These  simplifications  allow  us  to  calculate  the  array 
performance  without  specifying  the  carrier  frequencies  or  the  symbol 
rate. 

The  covariance  matrices  for  the  interference  and  the  noise  signals 
are  identical  to  those  given  by  (3.7)  and  (3.8)  for  BPSK  with  uii=w0. 

The  total  covariance  matrix  *,  is  determined  by  summing  $4*  and  $n. 
The  array  weights  and  output  power  levels  are  calculated  using 
Equations  (2.5)  and  (2.28)  through  (2.34). 


C.  FSK  DETECTOR  PERFORMANCE  WITH  CM  INTERFERENCE  AND  AWGN 

In  this  section  we  determine  the  probability  of  error  for  the  FSK 
detector  of  Figure  6.2  when  the  desired  signal  is  corrupted  by  both 
white  Gaussian  nose  and  CW  interference.  A  similar  problem  is  addressed 
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in  [26]  and  the  results  of  this  analysis  are  summarized  here. 

We  assume  that  the  composite  signal  at  the  detector  input  is  given 

by 


s(t)  =  cos(o)i<t  +  Yk)  +  /2Pi  cos  (cu0t  +  Yi)  +  n(t)  (6.9) 

where  <  =  0  or  1  depending  on  the  transmitted  symbol,  y^  is  the  phase  of 
the  desired  signal  at  the  detector  input,  y^  is  the  interference  phase 
(where  it  is  assumed  that  y^,  Yq  and  y^  are  i.i.d  ~  U[0,2ir]),  and  n(t) 
is  a  white  Gaussian  noise  signal  with  two-sided  PSD  n/2. 

We  assume  that  the  detector  BPFs  are  ideal,  each  with  bandwidths 
1/T.  Therefore  the  desired  and  interfering  signals  will  nearly  reach 
their  steady  state  values  at  the  end  of  the  symbol  interval.  We  assume 
that  the  passbands  of  the  two  filters  are  not  overlapping  and  we  neglect 
any  intersymbol  interference  or  crosstalk.  The  noise  produces 
(independent)  zero-mean  Gaussian  random  processes  with  variance  n/T  at 
the  output  of  each  of  the  BPFs. 

At  each  BPF  output  we  either  have  Gaussian  thermal  noise  alone  or 
the  sum  of  a  CW  signal  and  the  Gaussian  noise  signal.  The  probability 
density  function  of  the  amplitude  of  the  filter  output  has  a  Rayleigh 
distribution  in  the  former  case  or  a  Rician  [27]  distribution  in  the 
latter  case.  We  denote  the  output  of  the  "0"  and  "1"  envelope  detectors 
at  the  end  of  the  given  symbol  interval  by  Zg  and  respectively. 

We  first  assume  that  a  "0"  symbol  is  transmitted.  In  this  case  the 
"0"  filter  has  the  sum  of  two  sinusoids  (the  desired  signal  and  the 
interference)  plus  the  thermal  noise  at  its  input.  The  "1"  filter  sees 
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only  the  noise  signal  at  its  input.  We  denote  by  y  the  phase  difference 
between  the  two  sinusoids  in  the  "0"  filter  (i.e.,  y  =  y0  -y-O-  Tho 
(Rician)  PDFs  of  the  amplitudes  of  the  filter  outputs  conditioned  on  y 
are  then  given  by. 


,  ,  „  z0  A00z0 
f(z0lY,"0")  =  IQ  exp 


z0  +  A00 


(6.10) 


for  the  "0"  filter  and. 


f (Z1jy,"0") 


Z1  A10z) 

(n/T)  !0  (n/T)  exp 


(  Z)  *  Aio^ 

V  2"/T  J 


(6.11) 


for  the  "1"  filter,  where  Ig  is  the  modified  Bessel  function  of  the 
first  kind  of  order  zero  and  A-jj  is  the  amplitude  of  the  sinusoidal 
signal  at  the  input  to  the  i  filter  when  symbol  j  is  transmitted. 
Therefore,  Ago  and  A10  are  given  by, 

A00  =V2Pd  +  4VV i  cos^  +  2Pi  » 

and, 

A10  =  °*  (6.13) 

Since  no  sinusoidal  signal  is  present  at  the  input  of  the  "1" 
filter  when  a  "0"  is  transmitted  (i.e.  since  A10  =  0)  then  the  Rician 
distribution  given  in  (6.11)  is  actually  a  Rayleigh  distribution  [271, 
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f<ziiY>"0"’ ■'^7nexp  'Wm 


(6.14) 


We  now  assume  that  a  "1"  symbol  is  transmitted.  In  this  case,  each 
filter  has  both  a  CW  signal  and  a  Gaussian  noise  signal  at  its  input. 

The  (Rician)  PDFs  for  the  amplitudes  of  the  filter  outputs  are  given 


,  ,  x  z0  /Aoizo\ 

f  (Z0|y,’  °‘  )  =  (n/T)  xn  \(n/T)/ 


for  the  zero  filter  and. 


zn  +  adA 

2(n/T)  )  ’ 


(6.15) 


f (Z- Iy.-I") 


(n/T) 


/Anzn\ 
V  (n/T )J 


Z1  +  A11 
2(n/T)  , 


(6.16) 


where. 


A01  =  ^  2pi  > 


(6.17) 


AH  =  /2Pd  . 


(6.18) 


For  equally  likely  symbols,  the  probability  of  a  symbol  error 
conditioned  on  y  is  given  by, 


|  P[Z1  >  Z0|y,m0"  sent]  +  |  P[ZQ  >  Z^Y.'T’sent] 


(6.19) 


These  probabilities  can  be  directly  evaluated  [26]  to  give  us  the 


following  expression  for  the  conditional  symbol  error  probability. 
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P(e|y)  =  T 


WT  ’  /2n7T 


1  f A10A00 

2  l0  V  2n/T 


A2  +  A2 
*10  +  *00 


J  2n/T  ’  /2n/T 


A11A01 


-7( 


A2  +  A2 
*11  +  *01 


(6.20) 


where  Q(a,b)  is  Marcum's  Q-function  [24].  Since  A10=0  we  can  use  the 
relationships , 


0(0, b)  .  e-^2)2  , 


(6.21) 


ln(0)  =  1 


(6.22) 


to  simplify  (6.20).  The  resulting  expression  for  error  probability 
conditioned  on  y  is. 


11  f  1  Ao°  \  /  Aoi  A11 

P(e|y)  2  2  exp  l  "  2  /2n?T  /  +  Q  \  71 7/f  ’  Tln/T 


A01A11 


.2  A  *2 
1  A01  +  A11 
”  ?  (  ?n/T  ) 


(6.23) 
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We  integrate  P  (e | y )  over  the  (uniformly  distributed)  phase  angle  y 


to  obtain  the  symbol  error  probability. 


1  , 2* 

p<e>  / 


"  2  T0 


1  n00 


N  ( 


'01 


A01A11 
/2n/T 


exp 


A2  +  A2 
A01  +  A11 

2n/T 


'll 


2  /2^7T /  +  Q  V  /^/T  ’  /2n7T 


dy  < 


(6.24) 


This  is  the  general  expression  for  the  symbol  error  probability  for  the 
FSK  detector  with  a  single  CW  interferer  exactly  on  one  of  the 
signalling  frequencies. 

We  note  that  if  P.  =0  then  A.,  =  0  and  An_  =  A,  =  2/pT.  In  this 
i  ui  oo  li  □ 

case  we  can  use  (6.21)  and  (6.22)  to  further  simplify  (6.24).  In  this 
case  P(e)  becomes, 


P(e)  =-e 


(1/2)  (E.)/n 


(6.26) 

which  is  the  well  known  result  for  non-coherent  FSK  with  AWGN  but  no 
interference  [7]. 

The  general  expression  for  P(e)  given  by  (6.22)  was  evaluated 
numerically  to  determine  the  detector  performance.  Figures  6.3  and  6.4 
show  the  results  from  calculations  of  P(e)  for  several  SNR  and  INR 
values.  Simpson's  rule  was  used  to  evaluate  the  integral  in  (6.22) 
where  the  iterative  method  described  in  [24]  was  used  to  evaluate 
Marcum’s  Q  function.  Either  a  polynomial  or  an  asymptotic  approximation 
was  used  to  estimate  I0.  We  note  that  the  performance  of  the  binary 
FSK  detector  with  CW  interference  is  worse  than  that  of  either  of  the 
binary  PSK  detectors  discussed  in  chapters  III  and  V. 
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LOG  IP  JEM 


P  (E)  VS.  SNR 


Figure  6.3.  Performance  of  ideal  FSK  detector  with  CW  interference. 
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PIE)  VS.  INR 


Figure  6.4.  Performance  of  ideal  FSK  detector  with  CW  interference. 
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D.  PERFORMANCE  OF  THE  COMBINED  LMS  ARRAY/FSK  DETECTOR 


In  this  section  we  calculate  the  bit-error-probability  for  the  FSK 
system  that  includes  the  LMS  array.  As  noted  at  the  end  of  section  B  of 
this  chapter,  the  covariance  matrix  for  FSK  signalling  given  by  (6.8)  is 
identical  to  that  of  BPSK  (after  the  simplifying  assumptions  were  made 
in  both  cases).  Therefore,  the  adaptive  array  calculations  for  the  FSK 
and  BPSK  cases  are  the  same.  Thus  to  calculate  the  system  performance 
we  need  only  to  calculate  P(e)  for  the  FSK  detector  when  the  array 
output  signals  appear  at  the  detector  input. 

Figures  6.5  through  6.10  show  the  system  performance  for  0-j  values 
of  10,  20  and  80°.  The  array  bandwidth  factor  k  is  set  equal  to  5  in 
each  of  the  cases  shown.  Qualitatively,  the  FSK  curves  are  very  similar 
to  those  calculated  for  the  other  modulation  methods.  The  FSK  error 
probabilities  are  slightly  worse  than  either  of  the  binary  PSK  cases  and 
roughly  equal  to  those  to  QPSK  (for  the  same  0j  and  k  values).  Actually 
FSK  is  slightly  better  than  QPSK  in  regions  near  the  "humps"  in  the  P(e) 
vs.  INR  curves,  but  worse  in  the  other  regions.  Results  for  FSK  cases 
with  k=10  were  also  calculated  but  are  not  presented  here  since  they 
display  no  unexpected  characteristics. 
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LOG  IP  (EH 


P  (E)  VS.  SNR 


SNR  (08) 


Figure  6.5.  FSK  P(e)  vs.  SNR  for  3-element  array  (9-j=10°,  k=5). 


145 


7“ - 1 - 1 - 1  i  r  i 

•M.  -to.  -10.  0.  10.  tO.  90. 

I NR  (DB> 


Figure  6.6.  FSK  P(e)  vs.  INR  for  3-element  array  (9-j=10°,  k=5). 
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LOG  IP  IE)) 


P (E)  VS.  INR 


Figure  6.8.  FSK  P(e)  vs.  INR  for  3-element  array  (0^=20°, 
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P  (E)  VS.  SNR 
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CHAPTER  VII 


PERFORMANCE  OF  A  NARROWBAND  ADAPTIVE  ARRAY 
WITH  BPSK  SIGNALLING  AND  CW  INTERFERENCE 

A.  INTRODUCTION 

Chapters  III  through  VI  examined  the  performance  of  several 
communication  systems  that  included  LMS  arrays.  We  found  that  P(e)  for 
these  systems  increases  as  the  array  bandwidth  increases.  In  this 
chapter  we  more  closely  examine  the  effects  of  bandwidth  on  system 
performance. 

In  the  previous  chapters  we  found  that  the  system  performance  was 
similar  for  BPSK,  QPSK,  DPSK,  and  FSK  desired  signal  modulation. 
Therefore,  in  this  and  the  following  chapters,  we  limit  our  discussion 
to  the  BPSK  system  with  the  understanding  that  similar  qualitative 
performance  would  be  achieved  for  each  of  the  other  modulation  methods. 

In  Chapters  III,  IV,  and  V  we  used  correlator-type  PSK  detectors. 
The  P(e)  calculations  that  we  made  assumed  that  the  detector  input 
bandwidth  was  wide  enough  so  that  the  desired  signal  was  passed  with 
negligibly  small  distortion.  Therefore,  this  model  did  not  allow  us  to 
calculate  P(e)  for  values  of  k  less  than  approximately  5. 

In  this  chapter  we  describe  modified  BPSK  signalling  and  detector 
models  that  allow  us  to  calculate  P(e)  for  k  values  as  small  as  1/2. 
These  models  are  based  upon  the  Nyquist  pulse  shaping  criterion  [16,28]. 


151 


This  criterion  states  that,  with  appropriate  transmit  and  receive 
filtering,  we  can  communicate  over  a  bandlimited  channel  with  no 
intersymbol  interference  ( IS  I ) .  Using  this  Nyquist  criterion  we  can 
limit  the  transmitted  signal  bandwidth  to  (1/T)  Hz.  Furthermore,  for  a 
fixed  input  SNR,  the  bandlimited  BPSK  system  can  achieve  a  P(e)  value 
equal  to  that  of  standard  (infinite  bandwidth)  BPSK.  Several  practical 
aspects  of  narrowband  signalling  are  discussed  in  [29]. 

In  the  next  section  we  examine  the  performance  of  the  narrowband 
BPSK  system  without  an  adaptive  array.  In  section  C  we  calculate  the 
system  performance  when  the  adaptive  array  is  added  to  reject  CW 
interference.  In  Section  C  we  also  examine  the  effects  of  the  array 
input  bandwidth  on  this  performance. 

B.  BPSK  SIGNALLING  OVER  A  NARROWBAND  CHANNEL 

In  this  section  we  describe  a  narrowband  BPSK  communication  system. 
In  this  system  the  transmitted  signal  bandwidth  is  limited  to  (1/T)  Hz. 
This  bandwidth  corresponds  to  a  bandwidth  factor  k  (as  defined  in  the 
previous  chapters),  of  1/2.  At  the  receiver,  a  matched  filter  processes 
the  waveform.  The  output  of  this  filter  is  periodically  sampled  and 
symbol  decisions  are  made  based  upon  the  signs  of  the  output  samples. 

By  judicious  choice  of  the  transmit  and  receive  filters  we  can  eliminate 
ISI  by  insuring  that  the  receiver  response  to  each  symbol  is  zero  at  the 
sampling  instants  corresponding  to  all  other  symbols. 

Figure  7.1  shows  the  narrowband  BPSK  system  model.  The  transmitted 
baseband  signal  f(t)  is  passed  through  the  transmit  filter  which  has  a 
transfer  function. 


152 


153 


Figure  7.1.  Narrowband  BPSK  communication  system 


1 


(7.1) 


HT(U,)  =  P(1r/T)^) 


sinc(^I) 


where. 


P("/T)(“)  =  1  „ 


-IT  /T  <  0)  <  TT /T 
]  oi  |  >  tt/T 


(7.2) 


To  understand  how  this  filter  works  we  examine  its  response  to  a 
single  input  pulse.  Figure  7.2(a)  shows  a  simple  input  pulse 
f Figure  7.2(b)  shows  the  filter  transfer  function  Ht(uj). 
The  filter  output  (in  the  time  domain)  is  given  by. 


g(t)  =  f(t)  *  hT(t)  (7.3) 

where  h^(t)  is  the  transmit  filter  impulse  response  and  *  denotes 
convolution.  Application  of  the  convolution  theorem  of  Fourier  analysis 
yields, 


G(oi)  =  F(oi)Ht(ui)  (7.4) 

where  G(oi)  and  F(oi)  are  the  Fourier  transforms  of  g(t)  and  f(t) 
respectively.  We  calculate  F(oi)  using  a  standard  Fourier  transform 
pair. 


F(ui)  =  T  sinc(uiT/2).  (7.5) 

We  evaluate  G(ui)  from  (7.1),  (7.4)  and  (7.5), 

G(w)  =  Tp^T^(ui).  (7.6) 
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-  v/Z  w/Z 

(b)  Transmit  filter  transfer  function 


at 


G(w) 

- . - T 


—ir/T  w/T 

(c)  Fourier  transform  of  filtered  pulse 


(d)  Filtered  baseband  pulse 


It 

2 

1 

“wd 

"d 

(e)  Fourier  transform  of  carrier  modulated  by  a  single  filtered  pulse 
Figure  7.2.  Transmitter  processing  of  baseband  pulse. 
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Figure  7.2(c)  shows  6(oj).  Application  of  another  standard  Fourier 
transform  pair  yields, 

g(t)  =  sinc(irt/T).  (7.7) 

This  waveform  is  shown  in  Figure  7.2(d). 

As  we  have  just  shown,  the  transmit  filter  converts  the  rectangular 
baseband  digital  pulses  into  a  series  of  sine  pulses.  Note  that  the 
transmitted  symbols  are  separated  by  T  seconds  and  sinc(irt/T)=0  for 
t=nT  (n=l,  2,  3,  ...).  Therefore,  the  single  symbol  represented  by  f(t) 
will  produce  no  ISI  for  any  preceding  or  following  symbols. 

g(t)  is  multiplied  by  a  carrier,  cos(ajC|t),  prior  to  transmission. 
The  transmitted  signal  is  therefore  given  by, 

g(t)cos(cj<jt)  =  sinc(irt/T)cos(£oc|t).  (7.8) 

Figure  7.2(e)  shows  the  Fourier  transform  of  the  transmitted  signal. 

As  the  signal  propagates  through  the  channel,  it  is  attenuated  and 
corrupted  by  an  AWGN  process  n(t)  with  two-sided  PSD  N0/2.  Later  we 
consider  the  effects  of  interference  added  as  the  signal  propagates 
through  the  channel . 

At  the  receiver  the  signal  is  given  by 
TT  t 

r(t)  =  sinc(~)cos(aJdt)  +  n(t)  (7.9) 

where  we  have  neglected  the  propagation  delay.  This  signal  is 
multiplied  by  cos(u)<jt)  and  the  resulting  signal  is  passed  through  the 
receive  filter  which  is  an  ideal  low  pass  filter  with  cutoff  frequency 
ir/T  radians  (1/2T  Hz).  This  filter  removes  the  double  frequency  term 
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and  limits  the  thermal  noise  power.  The  signal  at  the  receive  filter 
input  is  given  by 

Ad  irt 

r(t)cos(aJdt)  =  Y~  sinc(~)[l  +  cos(2(odt)]  +  n(t)cos(o>dt)  (7.10) 

The  output  from  this  filter  consists  of  a  desired  signal  term  Xj  and  a 
noise  term  Xn.  From  (7.10)  we  see  that 

^d  irt 

Xd=  — sinc(— )  (7.U) 

since  the  receive  filter  eliminates  the  double  frequency  term. 

The  noise  signal  at  the  input  to  the  filter  n(t)cos  (u^t)  has 
(two-sided)  PSD  N0/4  or  equivalently,  a  single  sided  PSD  of  N0/2.  Since 
the  receive  filter  has  a  bandwidth  of  1/2T  Hz  Xn  is  a  zero-mean  Gaussian 
random  variable  with  variance  (N0/2) ( 1/2T)=N0/4T. 

We  note  that  X^  reaches  its  maximum  value  at  t=0.  The  detector 
output  is  sampled  at  t=0  to  detect  the  symbol  represented  by  f(t).  The 
receiver  output  is  sampled  at  points  where  t=nT  (n=l,  2,  3,  ...)  to 
dettect  the  symbols  that  follow.  Figure  7.3  shows  the  filter  output 
(without  the  noise  signal)  for  the  transmitted  symbol  sequence  1011. 

From  this  figure  we  see  that,  unlike  the  conventional  BPSK  waveform,  the 
narrowband  BPSK  signal  does  not  have  a  constant  envelope.  The  varying 
amplitude  of  the  envelope  for  narrowband  BPSK  requires  that  the 
transmitter  must  have  a  peak-to-average  power  ratio  greater  than  unity. 

The  receiver  makes  symbol  decisions  based  upon  the  sign  of  the 
receive  filter  output  at  the  sampling  instants.  Positive  output  samples 
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produce  "1"  decisions  and  negative  output  samples  produce  "0"  decisions. 
The  symbol  error  probability  is  given  by 

P(e)  =  \  P(Xd  +  Xn  <  Opl"  sent)  +  \  P(Xd  +  Xp  >  OpO"  sent) 

"  -  '  !  (7.12) 

From  the  symmetry  of  signal  space  we  can  show  that  the  two 
conditional  probabilities  in  this  expression  are  equal.  Therefore  the 
expression  for  P(e)  becomes, 

P(e)  =|  P(Xd  +  Xn  <  0  pi"  sent) 

-Ad 

=  P(Xn<— )  (7.13) 

where  Xn  is  a  zero-mean  Gaussian  random  variable  with  variance  N0/4. 
Therefore 

^  )•  (7.14) 

Eb,  the  energy-per-bit  in  the  desired  signal  component  of  r(t),  is 
given  by  Ad2T/2  (this  is  easily  shown  by  application  of  Parseval's 
theorem  to  the  desired  signal  component  in  Equation  (7.10)).  Therefore, 
Equation  (7.14)  can  be  rewritten, 

(7.15) 

This  expression  for  P(e)  is  the  same  as  the  well  known  result  for  ideal 
detection  of  an  unfiltered  (and  thus  infinite  bandwidth)  8PSK  signal  in 
the  presence  of  AWGN.  Equation  (7.15)  implies  that,  with  a  bandlimited 
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a)  Individual  waveforms  representing  the  sequence  1011 


(b)  Composite  waveform  for  the  sequence  1011 
3.  Baseband  signal  for  simple  binary  sequence. 
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system,  we  can  achieve  performance  equal  to  that  of  an  infinite 
bandwidth  standard  BPSK  system. 

The  results  presented  above  were  derived  under  the  assumption  that 
AWGN  is  the  only  corrupting  signal  at  the  receiver  input.  We  now 
calculate  the  system  performance  when  r(t)  also  includes  a  random-phase 
CW  interfering  signal  with  amplitude  Aj .  We  assume  that  the 
interference  is  at  a  frequency  where  it  is  not  rejected  by  the  receive 
filter.  In  [18]  Rosenbaum  shows  that,  in  this  case,  P(e)  is  given  by. 


P(e) 


I 


0 


erfc 


n/2 E7 

+  V  "NS-  C0StTrei ) 


(7.16) 


where  E.  and  y  ,  (which  can  be  considered  a  dummy  variable  of 
integration)  are  as  defined  in  Chapter  III.  The  derivation  of  this 
equation  is  similar  to  that  used  to  obtain  Equation  (3.44)  in  Chapter 
III. 

We  note  that  Equation  (7.15)  is  almost  identical  to  Equation 
(3.44),  our  general  expression  for  P(e)  for  BPSK.  The  only  difference 
is  the  factor  of  sinc(uiT/2)  in  the  interference  term  Equation  (3.44). 
This  term  is  caused  by  the  frequency  response  of  the  integrator  in  the 
integrate-and-dump  detector  used  in  Chapter  III.  The  frequency  response 
of  the  detector  that  we  use  in  this  chapter  is  constant  over  the 
detector  passband. 

The  preceding  analysis  suggests  that  we  can  transmit  (1/T)  bits  per 
second  over  an  RF  channel  with  a  bandwidth  (1/T)  Hz  with  the  same 
performance  as  we  would  have  with  standard  BPSK.  However,  to  achieve 


160 


these  results,  we  must  be  able  to  realize  the  transmit  and  receive 
filters  in  the  narrowband  model.  The  noncausal  nature  of  the  impulse 
response  of  these  filters  is  evidence  that  the  filters  are  not 
realizable. 

Good  approximations  to  ideal  Nyquist  pulse  shaping  filters  do 
exist.  Most  notably,  members  of  the  raised-cosine  family  of  filters 
[16]  provide  performance  equal  to  that  achieved  with  the  ideal  low-pass 
transmit  and  receive  filters  used  our  narrowband  model.  Actually,  the 
ideal  low-pass  filter  is  one  member  of  the  family  of  raised-cosine 
filters.  All  raised-cosine  filters  have  the  same  equivalent  noise 
bandwidth  as  the  ideal  filters  and  produce  and  same  P(e)  with  no  ISI. 
These  filters  have  a  passband  between  0%  and  1D0%  larger  than  the  ideal 
low-pass  filters  assumed  our  model.  Approximations  for  many  of  these 
raised-cosine  filters  can  be  realized  [30].  Reference  [16]  provides  a 
thorough  analysis  of  system  performance  with  non-ideal  filters  (without 
interference). 

C.  BPSK  SYSTEM  PERFORMANCE  WITH  AN  LMS  ARRAY  AND  CM  INTERFERENCE 

Figure  7.4  shows  a  block  diagram  of  a  narrowband  receiving  system 
that  includes  an  LMS  adaptive  array.  We  assume  that  the  transmitter  and 
channel  are  the  same  as  those  shown  in  Figure  7.1.  We  also  assume  that, 
as  in  previous  chapters,  the  desired  and  interfering  signals  arrive  from 
angles  9^  and  0-j  respectively  (measured  from  the  broadside  direction). 
The  BPFs  at  the  element  inputs  are  assumed  to  be  ideal  with  bandwidth 
k(2/T). 
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In  order  to  calculate  the  adaptive  array  performance  we  must 
calculate  the  ACF  of  the  desired  signal  d(t).  d(t),  the  analytic 

desired  signal  during  the  nth  symbol  interval,  is  given  by, 

d(t)  =  bnAdsinc(irt/T)exp(jajdt)  (r-1  )T  <  t  <  nT  (7.17) 

X  ^  XL 

where  b  =  +1  if  the  nL  but  us  a  1  and  b  =  -1  if  the  nir  bit  is  a  0. 
n  n 

Figure  7.5  shows  Sd(u>),  the  PSD  of  cf(t).  Rd(x),  the  ACF  of  3(t), 
the  inverse  Fourier  transform  of  Sd(w).  R~(t)  is  given  by, 

Rd(T)  =  Ad  sinc  (~)  exp(ja>d-r)  (7.18) 

We  can  use  Equations  (2.20)  and  (7.18)  to  calculate  the  desired 
signal  covariance  matrix.  If  we  make  our  usual  simplifying  assumption 
that  Td«T,  then  Rd(-r)  is  approximated  by, 

Rd(0  -  Ad  exp(jojdT)  (7.19) 

With  this  approximation  the  desired  signal  covariance  matrix  is  the  same 
as  that  which  we  used  for  each  of  the  previously  examined  modulation 
methods  (see  for  example  Equation  (3.9)). 

The  interference  and  noise  models  are  the  same  as  in  Chapter  III. 
Therefore,  since  the  desired  signal,  interference  and  noise  covariance 
matrices  are  the  same  as  in  Chapter  III,  the  adaptive  array  calculations 
for  the  narrowband  model  are  identical  to  those  for  the  wideband  model 
of  Chapter  III. 

As  discussed  in  the  previous  section  if  the  interference  is  at  a 
frequency  which  falls  within  the  receive  filter  bandwidth,  then  P(e) 
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will  be  given  by  (7.16).  As  previously  noted,  this  expression  is  the 
same  as  Equation  (3.44)  with  AtoT  set  equal  to  zero.  Although  the  signal 
and  receiver  structures  assumed  in  Chapter  III  and  the  present  chapter 
differ  significantly,  the  expressions  that  describe  the  performance  of 
both  the  array  and  detector  are  nearly  identical  for  the  two  cases.  In 
fact,  for  Ao)T=0,  the  mathematical  formulations  for  the  wideband  and 
narrowband  models  are  identical.  However,  the  range  of  k  values  over 
which  each  of  the  models  is  appropriate  is  not  the  same.  The  wideband 
model  of  Chapter  III  is  only  valid  for  k  values  greater  than 
approximately  5.  For  the  narrowband  model  we  can  make  k  as  small  as 
1/2.  For  k  values  greater  than  5,  the  narrowband  and  wideband  systems 
(with  AojT=0)  perform  identically. 

Figures  7.6  through  7.11  show  the  (narrowband)  system  performance 
for  k=l/2.  For  this  k  value,  there  are  no  humps  in  the  P(e)  vs.  INR 
curves.  Therefore  we  might  be  tempted  to  state  that  the  bandwidth 
corresponding  to  k=l/2  is  optimum.  However,  upon  a  careful  comparison 
of  Figures  7.6,  3.9,  and  3.15,  we  see  that  P(e)  for  k=l/2  is  actually 
higher  than  it  is  for  k=5  or  k=10  at  the  left  hand  sides  of  the  curves. 
Thus  we  make  the  unexoected  observation  that,  for  low  INR  values,  the 
system  performance  improves  as  the  noise  power  at  the  array  input 
increases. 

For  easy  comparison,  curves  showing  the  performance  of  both  the 
wideband  (k=10)  and  narrowband  (k=l/2)  systems  are  plotted  in  Figure 
7.12.  This  figure  shows  P(e)  vs.  INR  curves  for  0,=1O  degrees  and  an 
SNR  of  12  dB.  For  INR  values  above  approximately  10  dB  the  narrowband 
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system  performs  better  than  the  wideband  system.  However,  for  INR 
values  below  10  d8  the  wideband  system  performs  better. 

There  is  a  simple  explanation  for  this  behavior.  As  the  INR 
increases,  the  array  begins  to  null  the  interference  when  the 
interference  and  noise  powers  at  the  array  input  are  approximately 
equal.  Therefore  the  null  begins  to  form  at  a  lower  INR  value  when 
k=l/2  than  it  does  when  k=10.  Consider  the  cases  shown  in  Figure  7.12 
for  an  INR  value  of  5  dB.  Figures  7.13(a)  and  (b)  show  the  adapted 
antenna  patterns  for  k=10  and  k=l/2  respectively.  For  k=10,  the  array 
has  not  formed  a  null  in  the  interference  direction  and  the  desired 
signal  is  very  close  to  a  pattern  maximum.  When  k=l/2,  a  null  is  formed 
on  the  interference.  Furthermore,  since  the  desired  and  interfering 
signals  are  spatially  close,  the  array  no  longer  keeps  the  desired 
signal  near  a  pattern  maximum.  Therefore,  the  SNR  at  the  array  output 
is  lower  with  k=l/2  than  with  k=10  for  the  case  shown.  The  increase  in 
P(e)  due  to  the  SNR  reducation  with  k=l/2  is  greater  than  that  which 
occurs  due  to  the  residual  interference  power  at  the  array  output  with 
k-10. 

Figure  7.14  shows  the  system  performance  of  SNR=12  dB,  0^=10°,  and 
k  values  of  1/2,  1,  and  2.  The  curves  shown  in  this  figure  show  that 
there  is  no  k  value  that  is  optimal  for  all  INR  values.  The  curve  for 
k=l  lies  significantly  below  that  for  k=l/2  for  INR  values  below 
approximately  12  dB.  For  larger  INR  values  the  P(e)  for  k=l  is  only 
very  slightly  greater  than  that  for  k=l/2.  For  k=2  we  see  even  better 
performance  below  12  dB  but  a  hump  is  beginning  to  form  in  the  curve 
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Figure  7.7.  P(e)  vs.  INR  for  3-element  narrowband  array  (9h=0°,  9i=10°,  ■ 
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vs.  INR  for  3-element  narrowband  array  (0(1=0 


P  (E)  VS.  INR 


12.  P(e)  vs.  INR  for  SNR=12  dB  for  two  different  k  values 
(9d=0°,  01=10°). 
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(a)  k  -  10 

Figure  7.13.  Adapted  antenna  patterns  for  two  k  values.  SNR=12  dB, 
9d=0°,  9l-10°. 
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(b)  k  =  1/2 

Figure  7.13.  Continued. 
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performance  below  12  dB  but  a  hump  is  beginning  to  form  in  the  curve 
above  12  dB.  P(e)  vs.  INR  curves  were  plotted  for  several  other  values 
of  SNR  and  0^  and  similar  results  were  observed.  In  each  case  the  hump 
in  the  curve  appears  for  k  values  greater  than  about  2. 

In  previous  chapters  we  examined  the  performance  of  the  various 
systems  when  the  interference  was  at  a  frequency  where  the  detector 
performance  was  not  degraded  (but  where  the  interference  was  within  the 
array  input  BPF  passbands).  For  the  narrowband  system,  if  |Au|  is 
greater  than  ir/T,  the  interference  will  be  completely  rejected  by  the 
detector.  The  performance  for  this  case  will  be  the  same  as  we  found 
for  the  wideband  system  when  the  interference  frequency  was  at  a  null  in 
the  detector  frequency  response.  Therefore,  the  results  shown  in 
Figures  3.20  and  3.21  depict  the  performance  of  both  the  wideband  and 
narrowband  BPSK  systems  when  the  interference  is  at  a  frequency  where 
the  detector  is  not  affected. 

In  this  chapter  we  have  examined  the  performance  of  a  bandlimited 
BPSK  system  with  an  adaptive  array.  We  found  that  the  humps  in  the  P(e) 
vs.  INR  curves  are  only  seen  for  k  values  greater  than  approximately  2. 
We  found  the  best  performance  for  k  values  between  1/2  and  2.  Since  the 
performance  for  k  values  between  1  and  2  is  never  much  worse  than  that 
for  k  =  1/2  (and  sometimes  better)  we  conclude  that  there  is  little 
reason  to  expend  greater  cost  or  effort  to  reduce  k  below  2. 
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CHAPTER  VIII 


PROTECTION  OF  A  BPSK  SYSTEM  FROM  WIDEBAND  INTERFERENCE 

A.  INTRODUCTION 

In  the  previous  chapters  we  considered  the  performance  of  digital 
communication  systems  in  the  presence  of  CW  interference.  In  this 
chapter  we  determine  the  protection  that  the  adaptive  array  can  offer 
to  a  communication  system  with  bandlimited  AWGN  interference.  We 
calculate  the  performance  of  the  narrowband  BPSK  system  described  in 
Chapter  VII  as  a  function  of  the  signal  powers,  arrival  angles,  and  the 
interference  and  array  input  bandwidths. 

We  also  examine  the  performance  of  a  modified  LMS  array  with  noise 
interference.  This  modified  array  uses  tapped  delay  line  transversal 
filters  to  implement  the  element  weights.  These  filters  allow  the  array 
to  steer  nulls  that  are  much  wider  in  bandwidth  than  those  possible  with 
conventional  LMS  array.*  We  shall  find  that  the  modified  array  often 
provides  performance  better  than  the  conventional  array  for  AWGN 
i nterference. 

Section  B  describes  the  desired  signal  and  interference  models. 
Section  C  describes  the  results  of  performance  calculations  and 
discusses  the  effects  of  the  various  system  parameters.  Section  D 
summarizes  the  results  and  compares  these  results  with  those  of  previous 
chapters. 

*By  conventional  we  mean  the  LMS  array  with  a  single  complex  weight 
behind  each  element. 
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B.  SYSTEM  MODEL 


In  this  chapter  we  again  use  the  BPSK  modulation  method  described 
in  the  previous  chapter.  The  same  desired  signal  and  noise  models  are 
used;  however,  we  assume  that  the  interfering  signal  at  each  element 
input  is  a  zero-mean  Gaussian  random  process  with  power  (Pi ) i n -  Ti.e  PSD 
of  this  interference  is  shown  in  Figure  8.1.  The  center  frequency  is 
and  the  bandwidth  of  the  interference  is  1  where  the  prime  indicates 
that  we  are  considering  the  signals  at  the  input  to  the  element  BPFs. 

We  assume  for  simplicity  that  as  this  assumption  has  little 

bearing  on  the  derived  results.  We  shall  show  that  varying  slightly 
about  has  little  effect  on  the  array  performance.  Varying  tu-j  about 


(Si  ),n(w) 
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Figure  8.1.  Interference  PSD  at  element  BPF  inputs. 
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will  never  hurt  the  detector  performance  since  the  detector  is  simply 
a  multiplier  (where  the  array  output  is  multiplied  by  cos^t))  followed 
by  an  ideal  low-pass  filter.  In  fact  when  u^  differs  greatly  from  u^ 
then  the  detector  filter  may  completely  reject  the  interference. 

We  define  the  relative  bandwidth  of  the  interference  as, 

8i'  =  (Aoo.'Vok  .  (8.1) 

Similarly,  we  define  the  desired  signal  relative  bandwidth  by, 

Bd  =  (Au(|)/ud  ’  (8*2) 

where  Aui^  is  the  bandwidth  of  the  desired  signal  and  odj  is  its  center 
frequency.  The  relative  bandwidth  of  each  input  BPF  is 

Bin  =  2kBd,  (8.3) 

where  k  is  the  bandwidth  factor  as  described  in  the  previous  chapters. 

If  Bn  is  larger  than  B^n  then  the  input  BPFs  will  reject  some  of  the 
interfering  signal.  Therefore  the  relative  bandwidth  of  the 
interference  at  the  BPF  output  is 

Bi  =  MIN(B.n,Bi’)  (8.4) 

• 

where  MIN(A,B)  denotes  the  smaller  of  the  two  quantities  A  and  P.  Note 
that  we  could  not  write  B..  in  this  simple  form  without  the  assumption 
that  oj.  =  ojj.  The  bandwidth  of  the  interference  at  the  BPF  outputs  is 
given  by  Auk  ,  where 

Aw,.  =  ukB..  •  (8.5) 
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The  power  in  the  interference  at  each  BPF  output  is 

°i2  =  (Pi)inMIN(Bin/Bi,’1)‘  (8*6) 

The  ACF  of  ?(T)  is  given  by 

Ri  ( t )  =  o-2  sinc(Au^-r/2)exp(jft^-r).  (8.7) 

Substitution  of  Aw-j  from  (8.5)  into  (8.7)  yields, 

R~(t)  =  0i2  sine  [(0. 5)Bioii x3  expfjui^).  (8.8) 

From  Equations  (2.14)  and  (2.16)  we  see  that  4^  =  uj-T^  ;  therefore,  we 
have 

R~(t)  =  c^2  sine  [(0.5)Bi  (<|>i/Ti  )t]  exp  [j (^/^  )t],  (8.9) 

Now  we  can  substitute  this  expression  for  R~(t)  into  (2.21)  to  determine 
the  interference  covariance  matrix  With  this  substitution  4k 
becomes 

r  ^  ,  ~j^i  .  r  1 

sinc[-pi<j>i]e  smc[-Bi4>i]e 

2  .  ri,  i  J$i  .  r  L  l  "j'H 

*,  =  a-  sinc[-pi4,i]e  i  sinc[ -^B.^.Je 

,  j  2<J>  -j  .1  ,  j  <|>  "j 

_  sinc[Bi<|)i]e  sinc^^  Je  1 

(8.10) 

From  this  expression  we  can  justify  our  previous  assertation  that 
the  array  performance  is  not  very  sensitive  to  slight  variations  in  0)^ 


181 


about  Wfi.  The  terms  in  the  interference  covariance  matrix  represent 
Equation  (8.9)  evaluated  at  t  =  ±nT-j  where  n=0,  1,  or  2.  The  only 
variable  in  (8.9)  that  is  dependent  on  is  the  argument  of  the 
exponent  of  expCJ(4>i/T-j  )tH.* 

Using  (2.16)  we  can  expand  the  exponential  terms  in  (8.10)  to 
obtain, 

exp[±jn<j>. ]  =  exp[±jn(u)i/ajd)irsin0.]  .  (8.11) 

From  this  expression  we  see  that  small  fluctuations  in  0)^  about 
produce  small  fluctuations  in  the  interelement  phase  shift.  This  is  the 
same  effect  that  is  observed  when  8,  varies  slightly.  Since  a  change  in 
0-j  of  a  few  degrees  does  not  dramatically  change  the  array  performance 
we  can  assert  that  a  small  change  in  will  not  either.  This  assertion 
was  verified  by  performance  calculations  in  which  wd,  “i ,  Bd  and  B-j  were 
chosen  independently.  Very  slight  variations  in  performance  were  noted 
as  the  ratio  (ood/^-j )  was  varied  from  approximately  0.9  to  1.1. 

We  assume  that  <&d  and  $n  are  the  same  as  in  Chapter  VII.  The  array 
weights  can  be  calculated,  as  usual,  using  (2.9),  (2.25)  and  (2.5). 

Both  the  desired  signal  and  interference  are  modified  by  the  array. 
However,  we  presently  show  that  the  array  adds  very  little  distortion  to 
the  desired  signal  while  it  changes  the  spectrum  of  the  interference 
considerably.  The  array  acts  as  a  transversal  filter  when  processing  the 
signals  at  its  input.  Such  a  filter  is  shown  in  Figure  8.2.  The  delay 


♦Note  that  the  product  B^  in  the  argument  of  the  sine  function  is 
independent  of  u>i .  This  is  easily  shown  from  (8.1),  (8.4)  and  (2.16). 
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Figure  8.2.  Transversal  filter  with  two  delay  elements  and  three 
weights. 

between  filter  taps,  which  is  labeled  T  in  this  figure,  is  the 
interelement  time  delay.  The  transfer  function  of  the  filter  for  a 
three-element  array  is  given  by, 

H(oj)  =  wi  +  w2exp(-jwT)  +  w3exp(-j(u2T).  (8.12) 

If  9d  *  9^  the  interelement  time  delays  are  different  for  the 
desired  and  interfering  signals.  Therefore  these  two  signals  are 
processed  by  different  transfer  functions.  In  general,  the  LMS  weights 
produce  a  transfer  function  for  the  desired  signal  that  is  very  close  to 
the  ideal  all  pass  filter  (i.e.  H(u>)  has  a  constant  magnitude  and  linear 
phase  shift  over  the  frequency  range  of  interest).  The  LMS  weights 
typically  produce  a  transfer  function  for  the  interference  which  has 
nulls  at  frequencies  where  the  interference  PSD  is  large. 
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As  an  example  of  the  spectral  filtering  properties  of  the  3-element 
array,  the  equivalent  filter  transfer  functions  were  calculated  for  a 
specific  case.  Figures  8.3  and  8.4  show  the  desired  signal  and 
interference  transfer  functions  when  e<j=0o,  9i=10°,  SNR=12  dB,  INR=20 
dB,  B<j=0.1,  and  B-j  =0. 1 .  From  these  figures  we  see  that  the  array  has  a 
negligible  effect  on  the  desired  signal  while  it  significantly  modifies 
the  interference  signal  spectrum. 

The  power  of  the  desired  signal  at  the  array  output  is  given  by 
Equation  (2.32)  where  (t )  is  given  by  (7.17).  The  desired  signal  at 
the  array  output  is  processed  by  the  detector  which  consists  of  a 
multiplier  (where  the  signal  is  multiplied  by  cos^t))  and  an  ideal 
low-pass  filter  which  rejects  all  but  the  baseband  signal  components.  A 
power  loss  of  one  half  occurs  during  both  the  multiplication  and 
filtering  processes  so  the  desired  signal  power  at  the  detector  output 
is  Pq/4.  Therefore  the  desired  signal  at  the  detector  output  at  the  nth 
bit  sampling  instant  is  given  by 

Xd(n)  =  bn  71^721  sinc[(ir/T)(t-nT)]  (8.13) 

where  bn=l  if  the  n^  bit  is  a  "1"  and  bn=-l  if  the  n^  bit  is  a  "0". 

The  noise  signal  at  the  array  output  has  a  two-sided  PSD  of  n/2 
where  n/2  is  given  by  Equation  (3.14).  At  the  multiplier  output  this 
signal  has  a  two-sided  PSD  of  n/4  (at  baseband)  or  equivalently,  a 
one-sided  PSD  of  n/2.  The  detector  low  pass  filter  has  a  bandwidth  of 
1/2T.  Therefore,  the  noise  produces  a  Gaussian  random  variable  (GRV) 
Xn(n)  at  the  filter  output  at  the  n^  sampling  instant  with  variance. 
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Var[Xn(n)]  =  (n/2)(l/2T) 


=  n/(4T).  (8.14) 

Calculation  of  the  interference  power  at  the  detector  input  is  not 
as  straightforward  as  calculation  of  the  noise  power  is  even  though  the 
interference  signal,  like  the  noise  signal,  is  a  Gaussian  random 
process.  Since  the  thermal  noise  signals  are  statistically  independent 
white  noise  processes,  their  weighted  sum  at  the  array  output  is  also  a 
white  noise  process.  However,  the  interference  signals  are  correlated 
from  element  to  element;  therefore,  the  spectrum  of  the  interference 
signal  at  the  array  output  will  be  modified  by  the  interference  transfer 
function  of  the  array  as  shown  in  Figure  8.3. 

Since  all  of  the  processing  done  to  the  signals  is  linear,  the 
random  variable  produced  by  the  interference  signal  at  the  detector 
output  is  Gaussian.  Since  only  the  interference  in  the  frequency  band 
bounded  by  ±  Lu^/2  affects  the  detector  we  need  only  to  calculate  the 
array  output  power  of  the  interference  signal  in  this  band  to  calculate 
the  detector  performance.  We  denote  this  power  by  (Pi)det*  Calculation 
of  the  actual  PSD  of  the  interference  signal  at  the  array  output  is  not 
necessary. 

To  calculate  the  output  interference  power  (Pi  )det  we  u*-“  (2.33) 
with  R-j(t)  set  equal  to  the  ACF  of  the  portion  of  the  interference 
signal  at  the  array  input  that  is  in  the  frequency  band  of  interest. 

This  ACF  is  given  by. 
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(R~)det(t)  ■  Co^NINCUBj/B,)]] 


s1ncC( l/2)MIN(Bd ,B1  )<u1xn  expfj^t)  (8.15) 

<V  *>» 

Substitution  of  (Rj)det  for  Ri  1n  (2.33)  yields  (P^det  at  the  array 
output.  When  this  signal  is  applied  to  the  detector,  the  multiplication 
and  filtering  processes  reduce  the  interference  power  by  a  factor  of 
four  (just  as  it  did  for  the  desired  and  noise  signals).  Therefore  the 
interference  produces  a  GRV  at  the  detector  output  at  the  nth  sampling 
instant  with  variance, 

Var[Xi(n)]  =  (P, )det/4  •  (8*16) 

The  probability  of  an  error  for  the  nth  bit  is, 

P(e)  =  0.4P[Xi(n)  +  Xn(n)  <  -  I  "1"  sent] 

+  0.5P[X-j  (n)  +  Xn(n)  >  /(Pd/2)  |  "0"  sent]  (8.17) 

The  two  probabil ities  in  this  expression  are  equal.  Therefore  P(e)  can 
be  written 

P(e)  =  P[X-j  (N)  +  Xn(n)  >  /(Pd/2)~  |  "0"  sent].  (8.18) 

Since  X-j(n)  and  Xn(n)  are  independent  GRVs  with  variance  (Pi  )det/4  and 
n/4T  respectively,  their  sum  is  a  GRV  with  variance  (Pi)det/4  + 
Therefore  P(e)  is  given  by 
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Figure  8.4.  Frequency  response  of  LMS  array  in  the  interference 

signal  direction  (9d=0°,  01=10°,  SNR  =  12  dB,  INR  =  20  dB, 
k  =  1/2,  Bd  =  B1  =  0.1). 
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P(e)  =  erfc[/(P^72y  /  /(Pi )det/4  +  n/4T)j. 


(8.19) 


Multiplication  of  the  numerator  and  denominator  in  the  argument  of  the 
error  function  by  /3T  yields, 

P(e)  «  erfcC/2Eb/(ri+n)'j  (8.20) 

where  Eb  =  PdT  and  Ei  =  (Pi)detT* 

Since  both  the  interference  and  noise  signal  produce  independent 
zero-mean  GRVs  at  the  detector  output,  P(e)  is  dependent  only  upon  the 
SINR  at  the  dettector  input  where  we  make  the  definition, 

SINRdet  =  Ed/(Ei+n).  (8.21) 

With  this  definition  (8.19)  becomes, 

P(e)  =  erfc[/ (2  SlNRdet)]  •  (8.22) 

Note  that  if  Ej=0  then  P(e)  =  erfc[/ (2?^/^ )']  the  familiar  result  for 
BPSK  with  AWGN. 

C.  RESULTS  OF  PERFORMANCE  CALCULATIONS 

The  derivations  of  the  previous  section  can  be  used  to  calculate 
the  performance  of  the  BPSK  system  with  a  3-element  LMS  array  and  noise 
interference.  In  this  section  the  results  of  such  calculations  are 
presented  for  several  signal  scenarios. 

Figure  8.5  shows  the  results  from  typical  performance  calculations. 
This  plot  shows  P(e)  as  a  function  of  SNR  for  several  noise  bandwidths 
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for  8(1=0°,  9-j=10°,  and  k=l/2.  In  each  case  the  signal  to  interference 
power  ratio  (SIR)  is  -30  dB.  The  lowest  curve  in  this  plot  corresponds 
to  the  case  where  the  relative  bandwidth,  B-j '  =0.  The  other  curves 
correspond  to  B-j '  values  of  0.01,  0.025,  0.05  and  0.075.  From  these 
curves  we  see  that  the  array  offers  less  protection  as  the  interference 
bandwidth  increases. 

P(e)  was  calculated  as  a  function  of  SNR  for  several  different 
values  of  0-j  and  SIR.  Figure  8.6  shows  a  compilation  of  the  results 
which  allows  easy  comparison.  The  scale  of  each  of  the  small  plots  in 
Figure  8.6  is  the  same  as  that  of  Figure  8.5.  From  these  figures  we  see 
that  the  system  has  the  most  difficulty  in  suppressing  strong  broadband 
interference  when  the  desired  signal  and  interfering  signals  are 
spatially  close.  The  effects  of  a  large  interference  bandwidth  are  most 
detrimental  at  very  low  SIR  values  (i.e.  when  the  interference  is 
significantly  stronger  than  the  desired  signal  at  the  array  input).  For 
an  SIR  value  of  0  dB,  the  system  performance  with  broadband  interference 
(for  B-i  ’  <  0.075)  is  not  significantly  worse  than  it  is  with  B-j'  =0. 

In  general,  the  resolution  properties  of  a  linear  array  are 
dependent  upon  the  signal  arrival  angles.  A  linear  array  can  steer  a 
null  with  infinite  bandwidth  n  the  broadside  direction*.  However,  it 
becomes  more  difficult  to  achieve  a  broadband  null  as  the  interference 
arrival  angle  varies  from  the  broadside  to  the  endfire  direction. 


*  For  example,  a  3-element  linear  array  with  a  weight  vector 
W=[-1,2,-1]T  j,as  a  nun  -jn  the  broadside  direction  for  all 
frequencies. 
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Performance  of  a  3-element  array  with  noise  interference 
for  various  CIR  values  (Pj=0o,  k=l/2,  B(j=0.1).  values 
and  scales  are  the  same  as  those  of  Figure  8.5. 


In  order  to  determine  the  effects  of  arrival  angle  of  a  broadband 
interference  signal,  we  examine  several  cases  for  which  the  desired  and 
interfering  signals  are  separated  by  20°.  A  mosaic  of  the  resulting 
P(e)  vs.  SNR  plots  for  these  cases  is  shown  in  Figure  8.7.  This  figure 
shows  the  system  performance  for  two  combinations  of  arrival  angles  and 
three  SIR  values.  The  left-hand  curves  in  Figure  8.7  show  the 
performance  when  0-j=O°  and  0d  =  20°.  When  0j=O°  the  array  steers  a 
wideband  null  in  the  interference  direction  and  good  performance  is 
noted  at  all  SIR  levels. 

The  right-hand  plots  in  Figure  8.7  show  the  performance  when  0^=90° 
and  0d=7O°.  In  these  cases  the  array  has  significant  difficulty  nulling 
the  interference,  especially  for  the  broadband  interference.  Although 
the  performance  appears  poor  in  these  latter  cases  it  is  important  to 
note  that  P(e)  with  no  adaptive  array  would  be  nearly  0.5  for  each  of 
the  cases  considered. 

From  the  results  presented  above,  we  see  that  the  interference  is 
most  effective  when  it  is  relatively  strong  and  broadband.  In  each  of 
the  cases  that  we  have  examined  so  far  k  =  1/2  and  Bj‘  <  Bd.  We  next 
compare  the  performance  for  four  different  cases  including  some  cases 
where  B-j '  >  Bd  and/or  k  >  1/2.  The  cases  to  be  examined  are: 

(A)  Bd  =  0.1,  0  <  B-j '  <  0.075,  k  =  1/2 

(B)  Bd  =  0.01,  0  <  Bi’  <  0.075,  k  =  1/2 

(C)  Bd  =  0.1,  0  <  Bj *  <  0.075,  k  =  10 

(D)  Bd  =  0.01,  0  <  B-j '  <  0.075,  k  =  10 
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In  each  case  e<j=0°,  8i=10°  and  SNR  =  12  dB.  Case  (A)  is  the  case 
considered  previous  in  this  chapter  (see  Figure  8.6).  Case  (B)  is 
similar  except  that  the  desired  signal  bandwidth  is  reduced  by  a  factor 
of  10.  The  array  bandwidth  is  also  reduced  by  this  factor  since  k  =  1/2 
in  both  cases  (A)  and  (B).  In  case  (B)  some  of  the  interference  signal 
is  blocked  by  the  array  input  BPFs  when  B-j  1  is  greater  than  0.01.  k  is 

equal  to  10  for  cases  (C)  and  (D).  The  other  variables  are  the  same  as 

in  cases  (A)  and  (B)  respectively. 

Figure  8.8  shows  the  results  for  case  (A).  For  this  case  the 
interference  is  passed  by  both  the  input  BPFs  and  the  detector  filter. 
From  Figure  8.8  we  see  that  the  interference  bandwidth  has  little  effect 
on  the  system  performance  for  INR  values  below  about  15  dB.  For  large 
INR  values  the  system  performance  quickly  degrades  as  the  interference 
bandwidth  increases.  These  are  exactly  the  observations  that  we  made 
for  this  case  after  examination  of  Figure  8.6. 

Figure  8.9  shows  the  results  for  case  (B).  The  upper  two  curves  in 
this  figure  show  the  performance  when  B^'*0  and  B-j '=0.01.  For  these 
cases  B-j '  <  B^  and  the  performance  is  the  same  as  it  was  in  case  (A)  for 

these  B-j '  values.  For  the  other  three  cases  shown  in  this  figure 

'  >  Bjj.  In  these  cases  the  array  input  BPFs  reject  some  of  the 
interference  and  the  performance  monotonically  increases  as  8j ' 
increases.  Therefore,  for  the  narrowband  array,  a  fixed-power 
interference  signal  is  most  disrupting  when  its  bandwidth  is  equal  to 
the  array  bandwidth. 
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Figure  8.10  shows  the  results  for  case  (C).  As  in  case  (A),  the 
interference  in  this  case  is  passed  by  both  the  input  BPFs  and  the 
detector  filter.  We  see  from  Figure  8.10  that  the  performance  becomes 
bandwidth  dependent  for  INR  values  above  approximately  25  dB.  At  high 
INR  values  the  curves  quickly  diverge  indicating  poor  performance  with 
strong  wideband  interference.  The  only  difference  between  the  cases 
shown  in  Figures  8.8  and  8.10  is  the  difference  in  the  value  of  k. 
Examination  of  these  figures  indicates  the  increase  in  k  from  1/2  to  10 
produces  two  effects.  First,  for  the  case  when  k=10,  the  characteristic 
hump  apperas  in  the  P(e)  vs.  INR  curves.  The  appearance  of  this  hump  is 
not  surprising  since  we  have  seen  it  for  each  of  the  other  modulation/ 
interference  combinations  that  we  have  examined  so  far.  The  second 
important  difference  to  be  noted  when  comparing  Figures  8.8  and  8.10  is 
the  different  INR  value  at  which  the  curves  begin  to  diverge.  This 
value  is  approximately  15  dB  for  the  case  when  kl =1/2  and  25  dR  when 
k=10*. 

Figure  8.11  shows  the  results  of  the  calculations  for  case  (D).  In 
this  case  the  interference  is  passed  by  the  input  BPFs  but  is  partially 
blocked  by  the  detector  filter  for  Bj '  >  0.01.  For  the  cases  shown  in 
the  upper  two  curves,  the  interference  falls  entirely  within  the 
passband  both  the  input  and  detector  filters.  The  performance  shown  in 
these  two  curves  is  identical  to  that  shown  in  the  corresponding  curves 


*  Note  that  the  noise  power  at  each  element  input  is  13  dB  larger  for 
case  (C)  than  it  was  for  case  (A).  Therefore,  the  curves  begin  to 
diverge  at  approximately  the  same  value  of  P^ /Pn  at  the  array  input 
for  the  two  cases. 
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Figure  8.10.  3-element  array  performance  -  case  (C). 
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in  Figure  8.10.  However,  for  the  larger  values  of  interference 
bandwidth,  a  significant  portion  of  the  interference  is  rejected  by  the 
detector  filter.  For  these  cases,  the  performance  with  a  fixed-power 
jammer  monotonically  improves  as  the  bandwidth  increases. 

In  each  of  the  cases  discussed  above  the  system  performance  became 
worse  as  the  interference  bandwidth  increased  until  the  point  was 
reached  where  the  interference  and  desired  signal  bandwidths  were  equal. 
This  behavior  was  noted  for  each  of  the  array  input  bandwidths  examined. 
We  found  that  the  performance  was  best  for  small  k  values.  Therefore, 
for  best  performance,  the  desired  signal  bandwidth  and  the  array  input 
bandwidths  should  each  be  as  small  as  possible.  Similarly,  from  a 
different  point  of  view,  the  most  effective  jamming  strategy  is  to 
spread  the  interference  power  over  as  wide  a  bandwidth  as  possible  while 
keepig  the  bandwidth  less  than  or  equal  to  the  desired  signal 
bandwidth. 

It  is  interesting  to  compare  the  results  obtained  in  this  chapter 
for  zero-bandwidth  noise  interference  with  those  presented  in  the 
previous  chapter  for  CW  interference.  We  might  at  first  expect  these 
two  cases  to  produce  identical  results  since  the  adaptive  array 
calculations  are  identical  in  the  two  cases.  However,  the  equations 
from  which  P(e)  is  calculated  for  each  case  are  significantly  different. 
When  the  results  are  compared  for  the  two  cases  we  note  that 
significantly  poorer  performance  is  obtained  for  the  zero-bandwidth 
Gaussian  noise  case.  In  the  CW  interference  case  we  assume  that  we  have 
a  sinusoidal  interference  signal  with  a  fixed  amplitude  and  random 
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phase.  In  the  noise  case  we  assume  that  we  have  a  Gaussian  random 
process  .for  which  we  know  the  variance.  Figure  8.12  shows  a  typical 
member  of  the  ensemble  of  noise  waveforms  under  the  assumption  that 
the  bandwidth  is  not  zero  but  is  very  small*. 

For  maximum  effectiveness  with  a  fixed  average  power,  a  narrowband 
jammer  will  be  more  effective  if  the  CW  output  signal  is  modulated  with 
a  Gaussian  random  process  so  as  to  produce  an  output  waveform  similar  to 
that  shown  in  Figure  8.12.  If  such  a  signal  is  to  be  truly  Gaussian, 
the  transmitter  must  be  capable  of  transmitting  at  peak  power  levels 
which  greatly  exceed  the  average  power  level.  Whether  or  not  a 
particular  transmitter  can  do  this  depends  on  its  hardware  design.  If 
the  maximum  and  peak  power  capabilities  of  a  particular  narrowband 
transmitter  are  equal  then  the  transmitter  might  as  well  send  a  CW 
signal  of  constant  amplitude  so  that  the  maximum  power  level  is  always 
being  transmitted. 

0.  PERFORMANCE  OF  A  TAPPED  DELAY  LINE  LMS  ARRAY  WITH  WIDEBAND 
INTERFERENCE 

In  this  section  we  examine  the  performance  of  a  modified  LMS  array 
with  broadband  Gaussian  interference.  This  array  has  transversal 
filters  (i.e.,  tapped  delay  lines)  behind  each  element  instead  of  a 

*  If  we  were  to  assume  that  the  Gaussian  process  truly  has  zero 
bandwidth  then  the  interference  would  be  a  constant  amplitude  CW 
signal.  Our  Gaussian  assumption  would  therefore  apply  to  the 
statistics  of  the  ensemble  of  possible  interference  waveforms  and  not 
the  time  statistics  of  one  ensemble  member.  Therefore  the  zero- 
bandwidth  Gaussian  interference  process  is  not  ergodic. 
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Figure  8.12.  Sample  interference  waveform  with  -  0. 

single  complex  weight.  The  use  of  tapped  delay  lines  as  the  weighting 
elements  in  the  LMS  array  was  first  suggested  by  Widrow  [1].  Compton 
and  Rodgers  studied  the  performance  of  several  different  structures  of 
tapped  delay  line  (TOL)  LMS  arrays  in  [31],  We  first  present  a  short 
description  of  the  TOL  array  and  why  it  is  effective  against  wideband 
interference.  We  then  present  the  results  from  performance  calculations 
which  illustrate  the  performance  improvement  offered  by  this  type  of 
processing. 

We  first  examine  why  TDL  adaptive  array  processing  is  more 
effective  than  conventional  LMS  processing.  We  consider  a  2-element 
array  as  an  example.  A  2-element  array  has  a  transfer  function  for  the 
interference  of. 


H.  (cu)  =  w1  +  w2exp[-jo)iTi  ] 


(8.23) 


where  and  are  the  element  weights  and  T\  is  the  interelement  time 
delay  given  by  (2.14).  If  the  array  is  to  steer  a  null  on  the 
interference  H^co)  should  be  zero.  Thus,  for  a  null,  we  require, 

wx  =  ^expC-joj.T.]  .  (8.24) 

Since  w^  and  w^  are  simply  complex  numbers,  (8.24)  can  only  be 
satisfied  at  a  single  frequency  uj..  Thus,  with  appropriate  weights,  the 
2-element  array  can  completely  cancel  a  CW  interfering  signal.  However, 
complete  cancellation  is  impossible  for  broadband  interference  since  the 

element  weights  are  frequency  independent. 

0 

Now  let  us  examine  the  transfer  function  of  the  2-tap  TDL  filter 
shown  in  Figure  8.13.  The  transfer  function  of  this  filter  is, 

H(oj)  =  Wj  +  w2exp[-ju)T]  .  (8.25) 

where  T  is  the  time  delay  produced  by  the  delay  line.  If  we  use  such  a 
filter  to  implement  each  of  the  element  weights  in  an  adaptive  array  we 
might  expect  to  improve  the  performance  against  broadband  interference. 
This  improvement  is  to  be  expected  since  the  TDL  element  weights  are  not 
frequency  independent.  Therefore  the  TDL  filters  can  (approximately) 
compensate  for  the  frequency  dependence  of  the  complex  exponential 
term  in  (8.24). 

Compton  [32]  studied  the  optimum  delay  line  length  and  number  of 
taps  for  suppression  of  a  single  wideband  jammer.  The  results  of  this 
study  indicated  that  the  length  of  the  delay  line  is  not  particularly 
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Figure  8.13.  2-tap  transversal  filter. 

critical  as  long  as  it  is  shorter  than  approximately  one  wavelength.  In 
a  real  array,  extremely  short  delay  lines  should  be  avoided  since  they 
may  require  a  very  large  dynamic  range  for  the  weights.  Compton  also 
shows  that  a  single  wideband  noise  jammer  can  effectively  be  suppressed 
with  just  two  taps  per  delay  line. 

We  now  consider  the  3-element  TDL  LMS  array  shown  in  Figure  8.14. 
The  signal  received  by  each  array  element  is  processed  by  a  2 -tap  TOL 
with  complex  weights  and  a  quarter  wavelength  delay  between  taps.  The 
weighted  signals  are  summed  to  produce  the  array  output.  Conventional 
LMS  feedback  loops  are  used  to  drive  each  element  weight.  The  signals 
at  the  tap  outputs  are  given  by, 
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Figure  8.14.  3-element  tapped  delay  line  LMS  array. 
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x2(t)  =  d(t)  +  ?(t)  +  nj(t) 

(8.26) 

x2(t)  =  3(t-Td)  +  T(t-Ti)  +  n2(t) 

(8.27) 

x3(t)  =  d(t-2Td)  +  T ( t-2T i )  +  K3(t) 

(8.28) 

x4(t)  =  x^t-T) 

(8.29) 

x5(t)  =  x2(t-T) 

(8.30) 

xg(t)  =  x3(t-T) 

(8.31) 

where  Td  and  TV  are  the  interelement  time  delays  for  the  desired  signal 
and  interference.  The  desired  signal  covariance  matrix  is  given  by. 


Rd<°> 

Ri(-Td) 

Rd<V 

Ra(°) 

Ri(2Td> 

r;<v 

r;<t) 

Rj(-Td*T) 

r;<vt> 

r;<t> 

Rd<2Td+T> 

RiJ(VT> 

Ri(-2Td)  Ri<-T> 

r;<-v  Ri'Td-T> 

RJ(0)  r;(2t„-t) 

R^-2Td+T>  R^°> 

R^(-Td+T>  W 

Rd<T)  Ri<2Td> 


R^(-Td-T)  Ri(-2Td-T> 

Rd<-T) 

Ri<-Td-T> 

Rd<Td-T> 

Rd<-T> 

Rd'-Td> 

Rd<-2Td> 

Rd(°) 

Rd<-Td> 

R^{0) 

(S.32) 
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Similarly, 


R~(0) 

rT(-2Ti> 

R^(-T) 
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R~(0) 
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(8.33) 
by, 

Rn('T)  0  0 

0  R~(-T)  0 

0  0  Rp  ( -T ) 

R~(0)  0  0 

n  ' 

0  R~(0)  n 

n 

0  0  R~(0) 

n 

(8.34) 

where  we  denote  the  noise  ACF  at  each  element  by  R~(t). 

n 


and  the  noise  covariance  matrix  is  given 

r;(o)  0  0 

0  R~(0)  0 

n v  ' 

o  o  r;«j) 

$i  = 

R~(T)  0  0 

0  R~(t)  0 

n '  ' 

o  o  r;(t) 
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We  assume  the  same  desired  and  interference  signal  models  used 
previously  in  this  chapter.  We  also  assume  that  the  element  BPFs  have 
the  same  bandwidth  as  the  desired  signal  (i.e.  the  bandwidth  factor 
k=0.5).  Furthermore,  for  simplicity,  we  shall  only  consider  the  case 
where  the  interference  bandwidth  is  less  than  or  equal  to  that  of  the 
desired  signal. 

The  calculations  of  the  noise  covariance  matrix  are  somewhat  more 
complicated  in  this  case  than  for  the  previous  cases  that  we  considered. 
In  each  of  these  previous  cases  the  noise  covariance  matrix  was  diagonal 
(i.e.,  $n  =  o2 I).  In  these  cases,  to  calculate  the  covariance  matrix, 
we  only  needed  to  know  the  total  noise  power  (a2)  at  each  element  input 
and  not  its  PSD  or  ACF.  In  the  present  case  we  must  also  include  the 
non-zero  off-diagonal  terms  in  the  noise  covariance  matrix. 

We  assume  that  white  Gaussian  noise  is  present  at  the  input  to  each 
element  input.  The  ideal  BPFs  at  the  element  inputs  limit  the  noise 
spectrum  to  a  rectangular  passband  as  shown  in  Figure  8.14.  We  assume 
that  the  total  noise  pased  by  each  filter  has  power  a2.  Therefore,  the 
noise  complex  ACF  is  given  by, 

Rn(t)  =  a2sinc(Au^T/2)exp(ju!jT).  (8.35) 

We  assume  that  the  array  reference  signal  is  of  the  form, 

r(t)  =  3(t)//(Pd)1n  (8.36) 

where  (P^)^  is  the  desired  signal  power  at  each  element  input.  The 
reference  correlation  vector  S  Is  then  given  by. 
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Rj(0) 

Rd<Td> 

«?CT) 

Rd(Td+T) 

Rd(2Td+T5 


(8.37) 


From  the  equations  given  above  we  can  calculate  the  array 
performance.  The  array  weights  are  calculated  first  using 


-1. 


w  -  [*d  +  $1  +  *nrs. 


(8.38) 


The  power  levels  of  the  desired  signal,  interference,  and  the 
thermal  noise  at  the  array  output  are  then  calculated  using  (2.32) 
through  (2.34). 

In  Chapter  II  we  defined  the  SINR  at  the  array  output  as 


SINR=p+p.  (8.39) 
i  n 

Earlier  in  this  chapter  we  found  that  P(e)  of  the  ideal  detector  was 
dependent  only  upon  the  SINR  at  the  detector  input  where  we  defined 
SINRqgt  as 


E 


d 


SINR .  .  =  7T~i — 77 
det  'Ei  W+  n 

d  d1 »  ti  -  Ipi)detT  and  n  is 
the  detector  input.  We  defined  (pi )det 


(8.40) 

the  one-sided  real  noise  PSD  at 
as  the  power  in  the  portion  of 
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the  interference  signal  that  lies  within  the  band  of  frequencies  for 
which  the  detector  is  sensitive.  Since,  in  the  present  case,  we  assume 
that  the  interference  bandwidth  is  less  than  or  equal  to  the  desired 
signal  bandwidth  then  )det  =  Pi  *  T*1US  (®*^0)  becomes 


SINR 


det  P^T  +  n 


(8.41) 


or  equivalently 


SINR 


det 


+  n/T  * 


(8.42) 


Since  we  assumed  that  the  bandwidth  of  the  BPFs  at  the  element 
inputs  is  1/T  then  n/T  is  the  noise  power  Pn  at  the  detector  input 
(i.e.,  the  array  output).  Therefore 


SINR 


det 


Pi  +  Pn 


(8.43) 


A  comparison  of  (8.39)  and  (8.43)  indicates  that  SINR^ef  uP°n 
which  P(e)  is  a  simple  function  of,  is  equal  to  the  SINR  at  the  array 
output,  i.e.. 


SINRdet  =  SINR.  (8.44) 

Using  this  expression  in  (8.22)  we  have 

P(e)  =  erfc[(2)(SINR)].  (8.45) 

Thus  P(e)  for  the  case  at  hand  can  be  simply  calculated  from  the  SINR  at 
the  array  output. 
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We  now  use  the  derivations  described  above  to  calculate  the 


performance  of  the  TDL  LMS  array  with  BPSK  signalling  and  broadband 
noise  interference.  As  an  example  we  consider  one  of  the  cases  for 
which  we  found  particularly  poor  performance  with  the  standard  LMS 
array.  We  examine  the  system  performance  when  the  input  signal  scenario 
is  identical  to  that  for  the  case  shown  in  Figure  8.5.  For  this  case 
0d=°o,  9i=10°,  SIR  =  -30  dB,  and  B1  =  0,  0.01,  0.025,  0.05,  and  0.075. 
The  results  for  the  TDL  LMS  array  are  shown  in  Figure  8.15. 

A  comparison  of  Figures  8.5  and  8.15  indicates  that  adding  the  TDL 
processing  to  the  LMS  array  greatly  improves  performance.  In  fact,  the 
performance  for  all  of  the  non-zero  B^  values  for  the  TDL  array  is  as 
good  as  the  B-j=0  case.  Therefore,  the  TDL  array  has  completely 
eliminated  the  performance  degradation  caused  by  non -zero  bandwidth  of 
the  single  jammer.  This  performance  improvement  was  noted  for  a  wide 
variety  of  input  signal  scenarios. 

Rules  for  designing  TDL  LMS  arrays  did  not  exist  until  recently 
[23]  and  several  aspects  of  their  performance  are  yet  to  be  studied. 

For  example,  the  performance  of  the  TDL  LMS  array  in  the  presence  of 
multiple  broadband  jammers  has  yet  to  be  examined  in  detail.  However, 
the  results  presented  here  indicate  it  appears  that  continued  study  of 
the  TDL  LMS  array  would  be  worthwhile. 
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Figure  8.15.  Performance  of  the  3-element  TDL  LMS  array  for  the  same 
cases  shown  in  Figure  8.5. 
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CHAPTER  IX 


SUMMARY  AND  CONCLUSIONS 

In  this  dissertation  we  have  examined  the  performance  of  several 
communication  systems  that  use  IMS  adaptive  arrays  for  interference 
suppression.  P(e)  at  the  detector  output  was  used  as  the  measure  of 
system  performance. 

We  have  found  that  the  LMS  array  improves  the  communication  system 
performance  in  many  cases.  We  first  examined  the  performance  of  8PSK, 
QPSK,  OPSK  and  FSK  systems  with  CW  interference.  We  showed  that, 
qualitatively,  the  performances  of  these  systems  were  quite  similar.  We 
found  that  the  variations  in  performance  noted  among  the  different 
systems  with  the  same  input  signals  were  generally  consistent  with  the 
well  known  relative  performance  levels  of  their  respective  modulation 
methods.  We  found  that  no  particular  modulation  technique  provided 
unusually  good  or  poor  performance.  For  this  reason  we  used  BPSK 
modulation  in  the  latter  chapters  with  the  understanding  that  similar 
qualitative  performance  would  be  observed  with  any  of  the  other 
modulation  methods. 

We  found  that  the  system  performance  is  dependent  upon  several 
factors.  The  performance  is  best  when  there  is  a  large  spatial 
separation  between  the  desired  and  interfering  signals.  The  performance 
generally  becomes  worse  as  the  array  input  bandwidth  increases  beyond 
the  desired  signal  bandwidth.  We  defined  k  as  the  ratio  of  the  array 
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input  bandwidth  to  the  reciprocal  of  the  symbol  rate.  For  k  values 
greater  than  approximately  2,  humps  appear  in  the  P(e)  vs.  INR  curves. 
The  height  of  these  humps  becomes  larges  as  k  increases.  The  peaks  of 
the  humps  occur  at  INR  values  where  the  interference  power  is 
approximately  equal  to  the  thermal  noise  power  at  the  array  input. 

In  order  to  determine  the  optimum  system  bandwidth,  we  described  a 
minimal  bandwidth  BPSK  system  that  allows  the  BPSK  signal  to  occupy  a 
bandwidth  corresponding  to  k=l/2.  We  examined  the  performance  of  the 
narrowband  system  for  k  values  between  1/2  and  10.  In  general  we  found 
the  best  performance  for  small  k  values.  However,  for  weak 
interference,  we  found  that  the  very  narrowband  (i.e.,  k  =  1/2)  system 
performed  slightly  worse  than  systems  with  k=l  or  k=2.  We  found  that 
the  array,  with  k=l/2,  successfully  nulls  the  weak  interference;  but  the 
constraints  imposed  on  the  pattern  by  the  null  also  produce  a  slight 
decrease  in  the  desired  signal  power  at  the  array  output.  The  net 
effect  of  the  null,  in  this  case,  is  to  reduce  the  output  SNR  and  P(e). 
However,  the  differences  in  performance  noted  for  k  values  betwen  1/2 
and  2  are  smal 1 . 

We  calculated  the  performance  of  the  narrowband  BPSK  system  for  the 
case  where  the  interfering  signal  is  a  Gaussian  noise  signal  with 
non-zero  bandwidth.  We  found  that,  for  a  fixed  interference  power,  P(e) 
increases  as  the  interference  bandwidth  increases  until  the  interference 
and  desired  signal  bandwidths  are  equal.  Therefore,  the  interference  is 
most  effective  when  its  bandwidth  is  the  same  as  the  desired  signal.  We 
found  that  wideband  interference  was  especially  disruptive  when  the 
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input  interference  power  is  10  dB  or  more  larger  than  the  desired  signal 
power.  We  found  that  wideband  interference  is  most  effective  when  it 
arrives  from  a  direction  near  endfire  since  the  interelement  time  delay 
is  largest  in  this  direction.  We  suggested  that  a  nonlinear  element 
geometry  would  eliminate  this  susceptibility  to  endfire  broadband 
interference. 

We  next  calculated  the  effects  of  noise  interference  on  a  tapped 
delay  line  (TDL)  LMS  array.  In  this  array  transversal  filters  are  used 
to  implement  the  element  weights.  We  found  that  the  TDL  LMS  array  was 
not  sensitive  to  the  interference  bandwidth  for  any  of  the  signal  levels 
and  bandwidths  that  we  examined.  The  primary  disadvantage  of  the  TDL 
LMS  array  is  the  added  hardware  complexity. 

In  summary  we  found  that  the  LMS  array  can  offer  a  significant 
degree  of  interference  protection  to  digital  communication  systems.  The 
system  performance  is  not  unusually  dependent  on  the  modulation  method. 
For  best  results  we  found  that  the  desired  signal  and  array  input 
bandwidths  should  each  be  as  small  as  possible.  We  also  found  that  the 
system  performance  is  best  when  the  desired  signal  and  interference  were 
well  separated  spatially.  We  showed  that  wideband  interference  is  more 
effective  than  CW  interference  but  that  tapped  delay  line  processing 
greatly  reduces  the  system  susceptibility  to  wideband  interference. 
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